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ABSTRACT 

We investigate the well-known correlations between the dynamical mass-to-light ratio M/L 
and other global observables of elliptical (E) and lenticular (SO) galaxies. We construct two- 
integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/SO 
galaxies with SAURON integral-field stellar kinematics to about one effective (half-light) ra- 
dius R c . They have well-calibrated /-band Hubble Space Telescope WFPC2 and large-field 
ground-based photometry, accurate surface brightness fluctuation distances, and their observed 
kinematics is consistent with an axisymmetric intrinsic shape. All these factors result in an 
unprecedented accuracy in the M/L measurements. We find a tight correlation of the form 
(M/L) = (3.80 ± 0.14) x Oc/200 km g- 1 ) ' 84 ^ 07 between the M/L (in the I-band) mea- 
sured from the dynamical models and the luminosity-weighted second moment a c of the line- 
of-sight velocity-distribution within R e , The observed rms scatter in M/L for our sample is 
18%, while the inferred intrinsic scatter is ~ 13%. The (M/L)-a c relation can be included 
in the remarkable series of tight correlations between a c and other galaxy global observables. 
The comparison of the observed correlations with the predictions of the Fundamental Plane 
(FP), and with simple virial estimates, shows that the 'tilt' of the FP of early-type galaxies, 
describing the deviation of the FP from the virial relation, is almost exclusively due to a real 
M/L variation, while structural and orbital non-homology have a negligible effect. When the 
photometric parameters are determined in the 'classic' way, using growth curves, and the a c 
is measured in a large aperture, the virial mass appears to be a reliable estimator of the mass 
in the central regions of galaxies, and can be safely used where more 'expensive' models are 
not feasible (e.g. in high redshift studies). In this case the best-fitting virial relation has the 
form (M/L) v j r = (5.0 ±0.1) x R e al / (L G), in reasonable agreement with simple theoretical 
predictions. We find no difference between the M/ L of the galaxies in clusters and in the field. 
The comparison of the dynamical M/L with the (M/L) pop inferred from the analysis of the 
stellar population, indicates a median dark matter fraction in early-type galaxies of ~ 30% 
of the total mass inside one i? c , in broad agreement with previous studies, and it also shows 
that the stellar initial mass function varies little among different galaxies. Our results suggest a 
variation in M/L at constant (M / L) pop , which seems to be linked to the galaxy dynamics. We 
speculate that fast rotating galaxies have lower dark matter fractions than the slow rotating and 
generally more massive ones. If correct, this would suggest a connection between the galaxy 
assembly history and the dark matter halo structure. The tightness of our correlation provides 
some evidence against cuspy nuclear dark matter profiles in galaxies. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: kinematics and dynamics - galaxies: structure 
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1 INTRODUCTION 

Early-type galaxies display a large variety of morphologies and 
kinematics. Some of them show little rotation, significant kine- 
matic twists, and kinematically-decoupled components in their cen- 
tres. These systems are generally classified as ellipticals (E) from 
photometry alone, and are more common among the most mas- 
sive objects, which also tend to be redder, metal-rich and to ap- 
pear nearly round on the sky. The remaining early-type galax- 
ies show a well-defined rotation pattern, with a rotation axis gen- 
erally well aligned with the photometric minor axis. These ob- 
jects are often less massive, tend to be bluer and can appear 
very flat on the sky. These galaxies are gener ally clas sified either 
as E or as lenticulars (SO) from photometry iDavies et alJ ll983l: 
iFranx. Illingworth. & de Zeeuwll99ll : lKorm endv & Bender 1996). 

Despite the differences between individual objects, when con- 
sidering early-type galaxies as a class, it appears that they satisfy 
a number of regular relations between their global parameters. The 
first to be discovered was the relation L oc <j 4 between the lumi- 
nosity and t he velocity dispersion of the stars in elliptical galaxies 
iFaber & Jacksonlll976l) . These authors also already realized that, 
under simple assumptions, the observed relation implies a variation 
of the mass-to-light ratio (M/ L) with galaxy luminosity. 

It is now clear that the Faber- Jackson rel ation is th e projection 
of a thin plane, the Funda mental Plane (FP: iDressler et al1ll987l : 
iDiorgovski & Davii3ll987l) . which correlates three global observ- 
ables: the effective radius R c , the velocity dispersion a and the 
surface brightness J e at the effective radius. The relation is of the 
form R c oc a a lj: ■ If galaxies were homologous stellar systems 
in virial equilibrium, with constant M / L, one would expect a cor- 
relation of the form R c oc o 2 I~ x . The observed relation has in- 
stead the form R c oc cr 1,24 /^ 0,82 Ijorgensen. Franx. & Kiaergaard 
1996), with generally good agreement between different estimates 
and a weak dependence on the photometric band in the visual region 
jBernardi et all2003tlColless et all200ll and references therein). 

One way to explain the observed deviation, or 'tilt', of the FP 
from the virial predictions is to assume a power-law dependence of 
the M/L on the other galaxy structural parameters (mainly galaxy 
mass). The M/L variation could be due to differences in the stellar 
population or in the dark matter fraction in different galaxies. Other 
options are also possible however. The virial prediction for the FP 
is based on the assumption that galaxies form a set of homologous 
systems, both in the sense of having self-similar density distribu- 
tions and in terms of having the same orbital distribution. These 
spatial and dynamical non-homologies could then also explain the 
observed tilt of the FP. All these possible effects are known to oc- 
cur in practice, and they must contri bute, to som e degree , to t he 
observed tilt of the FP (e.g. Ivan Albada. Berlin. & StiavellH l 995 ) . 

Various attempts have been made to estimate the contri- 
bution of the possible effects, with sometimes contradictory 
results. A number of authors tried t o study the origin of the 
FP til t using approximate models. ICiotti. Lanzoni. & Renzinil 
1 1996) used theoretical arguments to show that non-homology, 
particularly in the surface-brightness profile of galaxies, could 
in principle play a major role in the tilt. Subsequent studies 
concluded that non-homology is indeed significant from the 
indirect argument that the expected change in M/L of the stellar 
population alone can not ex pl ain the observed M/L variation 
j^rag^^^^iMe^__^9^__|Pa^g^]oreoTOH^^^^Ca^^hc 

lFoAell^ormTarilMfcBTo\yi1 ~ 
1999). From the observed spatial non-homology, namely the 
inverse correlation of galaxy concentration and luminosity 



<Caon. Capaccioli. & D'Onofrid fl993l iGraham. Truiillo. & Caor] 
l2^)0*J) various works concluded that the contribution of the spatial 
non-homology of galaxies of different lumino sities is responsible 
for a s ignificant fraction of the observed FP tilt ( Prugni eT& Simierl 
1 19971: (Graham & Collessl Il997t iTruiillo. Burkert. & Belli l2004 . 
These authors constructed simple spherical isotropic models and 
used large samples of galaxies to estimate the variation in the 
virial M / L from the variation of the shape of galaxy profiles 
with luminosity. Although changes in the light-profile shape 
do not, on their own, sig nificantly affect the tilt of the FP (e.g. 
ITruiillo. Graham. & Caor3l200ll) they produce a variation in the 
velocity dispersion profile which may influence the measured 

M/L. 

Ivan der Marell ll99ll) . iMagorrian et alj 1 19981) and 
Gerha rd et alj 120011) constructed detailed dynamical models 
of smaller samples of galaxies, reproducing in detail the photome- 
try and long-slit spectroscopy observations. They also investigated 
the M/L and concluded that the variations observed in the models 
are consistent with the observed FP tilt. However, the fact that the 
spread in their measured M / L values was much larger than that in 
the FP, made the results inconclusive. 

In this work we revisit what was done in the latter studies, with 
some crucial differences: 

(i) The quality of our data allows a dramatic improvement in 
the accuracy of the dynami cal M/L d eterminations; for our galaxy 
sample SAURON jBacon et all 1200 ll Paper I) integral-field ob- 
servations of the stellar kinematics are available, together with 
HST/WFPC2 and ground-based MDM photometry in the I-band 
(which we reproduce in detail with our models including ellipticity 
variations), and reliable distances b ased on the s urface brightness 
fluctuations (SBF) determinations bv lTonrv et aljl200ll) : 

(ii) Our sampl e was extracted from the SAURON representative 
sample o Zeeuw et al. (2002, hereafter Paper II), which spans a 
wide range in velocity dispersion and luminosity, and includes both 
bright and low-luminosity E and SO galaxies; 

(iii) In contrast to previous studies the axisymmetric modelling 
technique we use can also be applied to very flattened galaxies, or 
objects containing multiple photometric components. 

All this allows us to carefully estimate the intrinsic scatter in the 
correlations between the dynamical M/L and other galaxy observ- 
ables, and makes it possible to set tight constraints on the measured 
slopes, for stringent comparisons with previous FP results and pre- 
dictions of galaxy formation theory. 

The paper is organized as follows. In Section 2 we present the 
sample selection and the photometric and kinematical data used. In 
Section 3 we describe the two-integral and three-integral dynamical 
models used to measure the M/L. In Section 4 we study the corre- 
lation of the M/L with global galaxy parameters. We discuss our 
results in Section 5 and summarize our conclusions in Section 6. 



2 SAMPLE AND DATA 
2.1 Selection 

The set of galaxies we use for this work was extracted from the 
SAURON sam ple of 48 E and SO galaxies (classification taken from 
Ide Vaucouleurs et alll99ll hereafter RC3), which is representative 
of nearby bright early-type galaxies (cz < 3000 km s _1 ; Mb < 
— 18 mag). As described in Paper II, it contains 24 galaxies in each 
of the E and SO subclasses, equally divided between 'cluster' and 
'field' objects (the former defined as belonging to the Virgo cluster, 
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the Coma I cloud, or the Leo I group, and the latter being galaxies 
outside these clusters), uniformly covering the plane of ellipticity, 
e, versus absolute blue magnitude, Mb- 

One of the requirements for the subsample defined in the 
present paper was the availabilit y of an accurate d istance deter- 
mined using the SBF technique bv lTonrv et alj J200lh . 38 objects in 
the SAURON sample satisfy this criterion. The second requirement 
was the availability of Hubble Space Telescope (HST) photometry 
obtained with the WFPC2 in the /-band (F814W). We selected this 
band because it is minimally affected by intrinsic dust effects and it 
is expected to trace well the bulk of the luminous mass in galaxies. 
35 galaxies satisfy this second criterion. The intersection between 
the two groups with SBF distances and J-band photometry leads 
to a sample of 29 galaxies. Out of these galaxies we additionally 
eliminated the objects showing strong evidence of bars within the 
SAURON field, from either the photometry (see Section l3~Tl or the 
kinematics. Specifically we excluded the five galaxies NGC 1023, 
NGC 2768, NGC 3384, NGC 3489 and NGC 4382. We did not at- 
tempt to eliminate all possibly barred galaxies from our sample, but 
at least we excluded the objects for which it makes little sense to 
deproject the observed surface density under the assumption of ax- 
isymmetry and to fit an axisymmetric model, which has to produce a 
bi-symmetric velocity field aligned with the photometry. We added 
the special galaxy M32, which is not part of the representative sam- 
ple, to the remaining sample of 24 galaxies. This last object allows 
us to sample the low-luminosity end of our correlations, which by 
construction cannot be accessed with galaxies from the main sur- 
vey. This leads to a final sample of 25 galaxies. 

2.2 Photometry 

Our photometric data consists of the HST/WFPC2/F814W im- 
ages already mentioned, together with ground-based photometry 
obtained, in the same F814W filter, with the 1.3m McGraw-Hill 
telescope at the MDM observatory on Kitt Peak. A relatively large 
field-of-view of 17.' 1 x 17f 1 was used for the MDM observations 
to provide good coverage of the outer parts of our galaxies and 
to allow for an accurate sky subtraction from the observed frames. 
The MDM images are part of a complete photometric survey of the 
SAURON galaxies, and they will be described in detail elsewhere. 

We adopted the WFPC2 images as reference for the photo- 
metric calibration and rescaled the MDM images to the same level. 
For this we measured logarithmically-sampled photometric profiles 
on the two images, using circular apertures, after masking bright 
stars or galaxies within the field. The sky level of the wide-field 
ground-based image can be determined easily and accurately by re- 
quiring the profile to tend asymptotically to a power-law at large 
radii. This leaves two free parameters: the scaling ratio between the 
ground-based image and the WFPC2 image, and the sky level of the 
WFPC2 image. The two parameters were fitted by minimizing the 
relative error between the two photometric profiles in the region of 
overlap, while excluding the innermost 3" to avoid seeing effects 
(top panel in Fig.0. After the fit we merged the two sky-subtracted 
and scaled profiles into a single one, which was used for the de- 
termination of the global photometric parameters (bottom panel in 
Fig.ffl- 

We measured the half-light radius R c and total galaxy lumi- 
nosity as in classic studies on the FP of e arly-type galaxies (e.g. 
iBurstein et alJll987l : lj0rgensen et al j fl*995). namely from a fit of a 
de Vaucouleurs R 1 ^ growth curve to the aperture photometry. The 
fit was generally restricted to radii R > 10" where the growth curve 
is a good description of the observations. We also experimented 
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Figure 1. Top Panel: match between the photometric profile of NGC 3379 
derived from the WFPC2/F814W (red solid line) and MDM images (dia- 
monds). The sky-level of the WFPC2 image and the scaling factor of the 
MDM image are the free parameters of the fit. The flattening in the central 
ground-based profile is due to saturation of the image inside R< 3". The 
fit was performed in the region between the two vertical dotted lines, where 
the WFPC2 and MDM profiles overlap, while excluding the innermost 3" 
to avoid PSF and saturation effects. Bottom Panel: fit of a de Vaucouleurs 
-R 1 / 4 growth curve (green solid line) to the combined WFPC2+MDM ob- 
served aperture photometry (diamonds). The two vertical dotted lines define 
the boundaries of the region that was included in the fit, which is limited to 
200" to reduce uncertainties due to the sky subtraction and contamination 
from nearby galaxies or foreground stars. The solid vertical line marks the 
location of R c . 




with alternatives: (i) fitting a lSersid fl9 68) R}/ n curve to the pro- 
file, (ii) fitting a R 1 ' n growth curve to the aperture photometry, (iii) 
fitting a R 1 ' n growth curve, with the value of n fixed by a previ- 
ous fit to the profile. We found that the values of R c obtained with 
different methods differ substantially in some cases. This is not too 
surprising as the determination is based on an extrapolation of the 
profile. We adopted the R 1 ^ growth curve fit because it allows a 
direct comparison of our results with previous works on the FP in 
Section l4~4l Another advantage is that with a growth curve approach 
the profile does not have to be well described by the assumed pa- 
rameterization, but one just tries to extrapolate the luminosity to 
infinite radius, starting from the outermost measured photometric 
points. This makes the measurement more robust and better repro- 
ducible than a profile-fitting approach. The measured values of the 
global photometric parameters are given in TableQ 
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Table 1. Galaxy sample and measured parameters. 



Galaxy Name 


R c 




It 


K T 




i 


(M/L) Joans 


(M/L) Srhw 


(M/L) pop 


(m - M) 


A B 


Fast-rotator ? 




(arcsec) 




(mag) 


(mag) 


(kms- 1 ) 


(°) 


(/-band) 


(/-band) 


(/-band) 


(mag) 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC 221 (M32) 


30 


0.92 


7.07 


5.09 


60 


50 


1.25 


1.41 


1.17 


24.49±0.08 


0.35 


yes 


NGC 524 


51 


0.61 


8.76 


7.16 


235 


19 


5.04 


4.99 


3.00 


31.84±0.20 


0.36 


yes 


NGC 821 


39 


0.62 


9.47 


7.90 


189 


90 


3.54 


3.08 


2.60 


31.85±0.17 


0.47 


yes 


NGC 2974 


24 


1.04 


9.43 


8.00 


233 


57 


4.47 


4.52 


2.34 


31.60±0.24 


0.23 


yes 


NGC 3156 


25 


0.80 


10.97 


9.55 


65 


67 


1.70 


1.58 


0.62 


31.69±0.14 


0.15 


yes 


NGC 3377 


38 


0.53 


8.98 


7.44 


138 


90 


2.32 


2.22 


1.75 


30.19±0.09 


0.15 


yes 


NGC 3379 (M105) 


42 


0.67 


8.03 


6.27 


201 


90 


3.14 


3.36 


3.08 


30.06±0.11 


0.10 


yes 


NGC 3414 


33 


0.60 


9.57 


7.98 


205 


90 


4.06 


4.26 


2.91 


31.95±0.33 


0.10 


no 


NGC 3608 


41 


0.49 


9.40 


8.10 


178 


90 


3.62 


3.71 


2.57 


31.74±0.14 


0.09 


no 


NGC 4150 


15 


1.39 


10.51 


8.99 


77 


52 


1.56 


1.30 


0.66 


30.63±0.24 


0.08 


yes 


NGC 4278 


32 


0.82 


8.83 


7.18 


231 


45 


4.54 


5.24 


3.05 


30.97±0.20 


0.12 


yes 


NGC 4374 (M84) 


71 


0.43 


7.69 


6.22 


278 


90 


4.16 


4.36 


3.08 


31.26±0.11 


0.17 


no 


NGC 4458 


27 


0.74 


10.68 


9.31 


85 


90 


2.33 


2.28 


2.27 


31.12±0.12 


0.10 


no 


NGC 4459 


38 


0.71 


8.91 


7.15 


168 


47 


2.76 


2.51 


1.86 


30.98±0.22 


0.20 


yes 


NGC 4473 


27 


0.92 


8.94 


7.16 


192 


73 


3.26 


2.91 


2.88 


30.92±0.13 


0.12 


yes 


NGC 4486 (M87) 


105 


0.29 


7.23 


5.81 


298 


90 


5.19 


6.10 


3.33 


30.97 ±0.16 


0.10 


no 


NGC 4526 


40 


0.66 


8.41 


6.47 


222 


79 


3.51 


3.35 


2.62 


31.08±0.20 


0.10 


yes 




14 


1 45 


10 40 


o.Vy 


1 10 


/ o 


2 81 


2 62 


1 44 




17 




NGC 4552 (M89) 


32 


0.63 


8.54 


6.73 


252 


90 


4.52 


4.74 


3.35 


30.87±0.14 


0.18 


no 


NGC 4621 (M59) 


46 


0.56 


8.41 


6.75 


211 


90 


3.25 


3.03 


3.12 


31.25±0.20 


0.14 


yes 


NGC 4660 


11 


1.83 


9.96 


8.21 


185 


70 


3.63 


3.63 


2.96 


30.48±0.19 


0.14 


yes 


NGC 5813 


52 


0.52 


9.12 


7.41 


230 


90 


4.19 


4.81 


2.97 


32.48±0.18 


0.25 


no 


NGC 5845 


4.6 


4.45 


11.10 


9.11 


239 


90 


3.17 


3.72 


2.96 


32.01±0.21 


0.23 


yes 


NGC 5846 


81 


0.29 


8.41 


6.93 


238 


90 


4.84 


5.30 


3.33 


31.92±0.20 


0.24 


no 


NGC 7457 


65 


0.39 


9.45 


8.19 


78 


64 


1.86 


1.78 


1.12 


30.55±0.21 


0.23 


yes 



Notes: (1) NGC number. (2) Effective (half-light) radius R c measured in the /-band from HST/WFPC2 + MDM images. Comparison with published values 
shows an rms scatter of 17%. (3) Ratio between the maximum radius /? max sampled by the kinematical observations and R c . We defined -R max = yj A/tt, 
w here A is the area on the sky sampled by the SAURON observations. (4) Total observed /-band galaxy magnitude. Comparison with the corresponding values 
by Tonrv et al. 12001) indicates an observational error of 13%. (5) Total observed /^T-band galaxy magnitude, from the 2MASS Extended Source Catalog 
I Jarrett et al. 2000). (6) Velocity dispersion derived with pPXF, by fitting a purely Gaussian LOSVD to the luminosity-weighted spectrum within R c . For the 
galaxies in which /J max / R c < 1 this value was determined on the whole SAURON field. Equation Q was used to derive the cr c value used in our correlations, 
adopting as R/R c the value in column 3. The statistical errors on these values are negligible, but we adopt an error of 5% to take systematics into account. 
(7) Inclination of the best-fitting two-integral Jeans model. We do not attach errors as they are fully model dependent. (8) M/L of the best-fitting two-integral 
Jeans model. (9) M/L of the best-fitting three-integral Schwarzschild model, computed at the inclination of column 7. Comparison with column 8 suggests 
an error of 6% in these values. (10) Stellar population (M / L) pop determined from measured line-strength values using single stellar population models. The 
median form al err or on this quantity is ~ 10%, but this value is strongly assum ption dependent. See Section l4"7l for details. (11) Galaxy distance mo dulus o f 
iTonrvetalJfeOOlh . adjusted to the Cepheid zeropoint of Freedman et al. 1 2001 ) by subtracting 0.06 mag (for a discussion see Section 3.3 o flMei et all20o3l . 
(12) 73-band galactic extinction of Schlegel, Finkbeiner, & Davis 1 1998) as given by the NED database. We adopt Aj = Ag /2.22 as in NED. (13) Galaxy 
classification according to the measured value of the specific stellar angular momentum within one R c , from the SAURON kinematics. A qualitative distinction 
between the two classes of galaxies can be seen on the kinematic maps of Paper III. More details will be given in Emsellem et al. (in preparation). 



2.3 Kinematics 

The kinematic measurements used here were presented in 
lEmsellem et alj 12004. hereafter Paper III), where further details 
can be found. In brief, the SAURON integral-field spectroscopic ob- 
servations, obtained at the 4.2-m William Herschel Telescope on 
La Palma, were reduced and merged with the XSauron software 
developed at CRAL (Paper I). They where spatially binned to a 
minimum signal-to-noise ratio S/N ~ 60 usin g the Voronoi two- 
dimensional (2D) binning algorithm by CaoDellari & Cooi il l2003l) 
and the stellar kinematics were subseque ntly extracted wi th the 
penali zed pixel-fitting (pPXF) method of ICappellari & Emselleml 
12004). This provided, for each Voronoi bin, the mean velocity V, 
the velocity dispersion a and four hi gher order Ga uss-Hermi te mo- 
ments of the velocity up to h^-ha Ivan der Marel & Franxlll993l: 
lGerhardf!99i . We estimated errors on the kinematics using 100 
Monte Ca rlo realizations an d applying the prescriptions of Sec- 
tion 3.4 of CaDDellari & Emsellem (2004). In addition we observed 



M32 with SAURON in August 2003 with the same configuration as 
for the other galaxies in Paper III, and obtained two pointings re- 
sulting in a field-of-view of about 40" x 60". 

For each galaxy in our sample we determined a e , the 
luminosity-weighted second moment of the line-of-sight velocity 
distribution (LOSVD) within the half-light radius R c . Compared to 
the central velocity dispersion <t c , which was sometimes used be- 
fore, this quantity has the advantage that it is only weakly dependent 
on the details of the aperture used. The observed a c is an approx- 
imation to the second velocity moment which appears in the virial 
equation iBinnev & Tremainell987l Section 4.3). It is proportional 
to vM (with M the galaxy mass), and so is weakly dependent on 
the details of the orbital distribution. These are the reasons why a 
similar quantity (estimated however from long-slit data) was also 
adopted by e.g. lGebhardt et all J2003) to study the correlation be- 
tween the mass of supermassive black holes (BHs) and the velocity 
dispersion of the host galaxy. 

The use of integral-field observations allows us to perform 
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the o e measurements in a more rigorous manner than was possi- 
ble with long slits. We measured o c from the data by coadding all 
luminosity-weighted spectra within _R C . The resulting single 'effec- 
tive spectrum' has extremely high S/N > 300 per spectral element, 
and is equivalent to what would have been observed with a single 
aperture of radius R c centred on the galaxy . We then used pPXF, 
with the same set of multiple templates from lVazdekisl 1 19991 here- 
after VZ99) as in Paper III, to fit a purely Gaussian LOSVD from 
that spectrum to determine o (this time the higher order Gauss- 
Hermite moments are set to zero and are not fitted). This proce- 
dure takes into account in a precise manner the effect of the galaxy 
rotation on o c , so that no extra correction is required. 

As we do not sample all galaxies out to one _R C we want to cor- 
rect for this effect and to estimate how much uncertainty this can in- 
troduce in our o a measurements. For this we show in Fig.|2|the pro- 
files of or, measured by coadding the SAURON spectra within cir- 
cular apertures of increasing radius, for all the 40 early-type galax- 
ies in Paper III which we sample to at least R c /2. In some cases the 
selected circular aperture is not fully sampled by our data, and we 
define the radius R = -J 'A/tt as that of a circular aperture with the 
same area A on the sky actually covered by the SAURON data inside 
that aperture. We fitted the profile of o e versus R for every galaxy 
with a power-law relation or oc R? and we determined the bi- 
weia ht {Hoaeli n. Mosteller. &Tukevl 19831) mean and standard de- 
viation for the exponent 7 = —0.066 ± 0.035 in the entire set. In 
addition we normalized the as values to the corresponding o c / 2 , 
defined as the dispersion at R c /2, and computed the biweight mean 
from all galaxies, at different fractions of R c . This mean or profile 
appears to be well described by the mean power-law relation fitted 
to the individual profiles. One can see that generally or decreases 
by less than 5% from R a /2 up to R c , but the galaxy to galaxy vari- 
ations are significant. In summary the aperture correction has the 
form: 



(or/o c ) = (R/R c 



-0.066±0.035 



(1) 



The logarithmic slope of our best-fitting relatio n is in reas onable 
agreement with the value —0.04 found by Jeirgensen etalj fl995t) 
and with the value -0.06 derived by iMehlert et alj J2003l) from 
long slit spectroscopy. 

The measured values of o c for our sample are given in TableQ 
together with the actual fraction of R c sampled by our kinematics. 
We decided not to restrict our o a measurements to a smaller aper- 
ture so as not to discard data for the galaxies that we sample to large 
radii. In what follows we correct the values of o c with the aperture 
correction of equation Q to determine our correlations and to gen- 
erate our plots. 



3 DYNAMICAL MODELLING 
3.1 MGE mass model 

We constructed photometric models for all the 29 candidate galax- 
ies (see beginning of Section l2~Tl plus M32 with HST photometry 
and SBF distances. For this we used the Multi-G aussian Ex pansion 
(MGE) parameterization by lEmsellem. Monnet. & Baconl ll994l) . 
which describes the observed surface brightness as a sum of Gaus- 
sians, and allows the photometry to be reproduced in detail, includ- 
ing ellipticity variations and strongly non-elliptical isophotes. The 
accurate MGE modeling of such a large sample of galaxies, each 
consisting of three separate images, was made feasible by th e use of 
the MGE fitting method and software by CaDDellari (2002), which 
was designed with this kind of large-scale application in mind. 
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Figure 2. Luminosity- weighted second moment of the line-of-sight velocity 
distribution within an aperture of radius R, normalized to its value at R, c / 2 
(see text details). This plot shows, with the thin coloured lines, all the 40 
E/S0 galaxies in the SAURON sample for which _R c /2 is not larger than the 
field-of-view. The open circles represent the biweight mean of all the curves 
at every given radius. The number of measurements sampled at every radius 
is not constant, moreover there are significant galaxy-to-galaxy variations 
in the profiles. The thick black line is a power-law relation ur oc R~ 066 
which has an exponent defined by the biweight mean of the individual values 
for every galaxy. The dashed lines indicate the standard deviation in the 
individual slopes. The galaxy with the largest increase of the ctr profile at 
large radii is NGC 4550 (see Paper III). 



Our MGE models were fitted simultaneously to three images: 
(i) the WFPC2/PC1 CCD (ii) the lower-resolution WFPC2 mosaic 
and (iii) a wide-field ground-based MDM image (Fig.[3J. In Sec- 
tion |^| we described how the ground-based and WFPC2 images 
were carefully sky-subtracted and matched. The MGE fits were 
done by keeping the position angle (PA) for all Gaussian fixed 
and taking PSF co nvolution into account (our F814W MGE PSF is 
given in Table 3 of CaoDellari et al J2002I) . The quality of the result- 
ing fits was inspected, together with the kinematics, to exclude the 5 
galaxies which could not be reasonably well fitted by a constant-PA 
photometric model (Section f2.lt . This reduced the original sample 
to 25 galaxies. 

As discussed in Caooellari (2002) we 'regularized' our solu- 
tions by requiring the axial ratio of the flattest Gaussian to be as 
big as possible, while still reproducing the observed photometry. 
In this way we avoided artificially constraining the possible incli- 
nations for which the models can be deprojected assuming they are 
axisymmetric. This also prevents introducing sharp variations in the 
intrinsic density of the MGE model, unless they are required to fit 
the surface brightness. This regularization of the models is needed 
because of the non-unique ness of the de projection of an axisym- 
metric density distribution (Rvbicki 1987). 

The resulting analytically deconvolved MGE model for 
each g alaxy was corrected for extinction following Schlegel et al. 
1 1998), as given by NED. It was converte d to a surface den sity in so- 
lar units using the WFPC2 calibration of Dolphin] feOOd) while as- 
suming an /-band absolute magnitude for the Sun of M t q = 4.08 
mag (all the absolute magnitudes of the Sun in this paper, in the 
I, R and A'-band, are taken from Table 2. 1 of Binnev & M errifiel3 



1998). The calibrated MGE models are given in Appendix B. The 
MGE model for NGC 2974 was taken from|lCrainovic et alj<2005l) . 
which uses precisely the same method and calibration. 
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Figure 3. Contours of the surface brightness of NGC 3379 (in 0.5 mag 
arcsec -2 steps) at three different scales. Top Panel: 35"x35" PCI WFPC2 
CCD. Middle Panel: the whole 160" X 160" WFPC2 mosaic. Bottom Panel: 
17' 1 X 17' 1 MDM image. North is up and east is to the left. The location on 
the MDM field of the L-shaped WFPC2 mosaic is also shown. The galaxy 
at the left of the frame is NGC 3384, which was masked during the MGE 
fit. Superposed on the three plots are the contours of the constant-PA MGE 
surface brightness model, convolved with the PSF of each observation. The 
MGE model was fitted simultaneously to all three images. 



3.2 Two-integral Jeans modeling 

We constructed axisymmetric Jeans models for the 25 galaxies in 
our sample in the following way. Under the assumptions of the 
MGE method, and for a given inclination i, the MGE surface den- 
sity c an be deprojected analytically l Mon net Bacon. & Emselleml 
1992) to obtain the intrinsic density p(R,z) in the galaxy merid- 
ional plane, still expressed as a sum of Gaussians. Although the 
deprojection is non-unique, the MGE deprojection represents a rea- 
sonable choice, which produces realistic intrinsic densities, that re- 
semble observed galaxies when projected from any inclination. 

For an axisymmetric model with constant M/L, and a stel- 
lar distribution function (DF) which depends only on two inte- 
grals of motion / = f(E,L z ), where E is the energy and 
L z the angular momentum with respect to the z-axis (the axis 
of symmet ry), the second velocity mom ents are uniquely de- 
fined (e.g. lLvnden-B"eill96l iHunteill 19771) and can be computed 
by solving the Poisson and Jeans equations (e.g. ISatohl fl98Cl : 
iBinnev. Davies. & Illingworthlll990|) . The second moment ji 2 pro- 
jected along the line-of-sight of the model is then a function of 
only two free parameters, i and M/L. We do not fit for the mass 
of a possible BH, as the spatial resolution of our data is not suffi- 
cient to constrain its value, but we set the ma ss of the centr al BH 
to that predicted by the M bh-Q" corr elation iFerrarese & Merritt 
l200rilGebhardt et all20 00) as given bv lTremaine et aljj2002l) . The 
inclusion of the BH is not critical in this work, but it has a non- 
negligible effect for the most massive galaxies. E.g., in the case of 
M87 the expected BH of ~ 3 x 1O 9 M prod uces a small but de- 
tectable eff ect on the o bserved a up to a radius ~ 5" (e.g. Fig. 13 
of lvan der Marelll994l) . To make our results virtually insensitive to 
the assumed BH masses we do not fit the innermost 7? < 2" from 
our data, where the BH may dominate the kinematics. 

Starting from the calibrated parameters of the fitted MGE 
model and using the MGE formalism, the model fi2,*, without a 
BH, can be computed easily and accurately usin g a single one- 
dimen sional integral via equations (61-63) 1 of Emselle m et alj 
ll994l) . The contribution /i2,. to the second m oments, due to 
the BH , is computed via equations (49,88,102) of Emselle m et alj 
1 19941) and also reduces to a single one-dimensional integral. Given 
the linearity of the Jeans equations with respect to the potential, the 
luminosity weighed second moments for an MGE model with a BH 
can be obtained as Vpfi 2 = v p fM 2ti , + fp/ttl,., where v v is the de- 
convolved MGE surface brightness. The u p fi 2 has to be convolved 
with the PSF and integrated over the pixels before comparison with 
the observables. For the 25 galaxies in our sample we computed the 
model predictions for fj, 2 , at different inclinations, for each Voronoi 
bin on the sky. At every inclination the best fitting M/L is obtained 
from the simple scaling relation {M/L) oc fi 2 as a linear least- 
squares fit: 



(M/L) 



d d 
d • m 



(2) 



where the vectors d and m have elements d n = fj,' 2 n /Afj.' 2 n and 
m n — /i2,n/A/i 2 . n , with fi' 2i n = VV„ + o'n the N measured val- 
ues and A^i 2 ,n me corresponding errors, while fi 2<n are the model 

predictions for (M/L) = 1. 

As previously noticed by other authors (e.g. Ivan der M arel 



Their equation (63) has a typographical error and should be replaced by 

^(l-e 2 T2) 



the following B ■ 
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Figure 4. Relative variation of M / L predicted from the virial equations as 
a function of inclination, for MGE models with different constant observed 
axial ratios q' going from an observed flat galaxy q' = 0.3 (green line) to an 
observed round object q 1 = 0.9 (red line). At each observed axial ratio the 
density is assumed to be rounder than an extreme value of q = 0.1. Note 
the small variation < 10% of the M/L estimate at different inclinations, 
for q' < 0.7. 



1991 ) the M/L derived from Jeans modeling is only weakly depen- 
dent on inclination, due to the fact that the increased flattening of a 
model at low inclination (i = 90° being edge-on) is compensated 
by a decrease of the observed velocities, due to projection effects. A 
quantitative explanation can be obtained by studying how the M/L 
determined from the tensor virial equations varies with inclination, 
for given observed surface brightness an d H2- This variation for an 
MGE model is given by equation (77) of Emsell em et alJll994l) . If 
the density is stratified on similar spheroids, the M/L variation is 
independent of the radial density profile, while only being a func- 
tion of the observed axial ratio of the galaxy and the assumed incli- 
nation. A plot of the variation of M/L with inclination, for differ- 
ent observed axial ratios, is presented in Fig.[4]and shows that one 
should expect errors smaller than < 10% even for significant errors 
in the inclination, if the observed galaxy has an observed axis ratio 
q' less than 0.7. When the galaxy is rounder than this, it is crucial 
to obtain a good value for the inclination, otherwise the measured 
M/L can be significantly in error. 

The best-fitting inclination provided by the two-integral mod- 
els needs not necessarily be correct, as there is no physical reason 
why the real galaxies should be well described by two-integral mod- 
els. Nonetheless we found that for the four galaxies in our sample 
which are not already constrained by the photometry to be nearly 
edge-on, and for which an independent estimate of the inclination 
is provided by the geometry of the dust or by the kinematics of the 
gas (NGC 524, NGC 2974, NGC 4150 and NGC 4459), the inclina- 
tion derived from the Jeans model agrees with that indicated by the 
dust or gas (see Appendix A). Moreover, the round non-rotating gi- 
ant elliptical galaxies (e.g. NGC 4552 and NGC 4486) are best fitted 
for an inclination of 90° , indicating that they are intrinsically close 
to spherical. This is in agreement with the fact that this whole class 
of massive galaxies, with flat nuclear surface-brightness profiles, 
which are often found at the centre of galaxy clusters, always ap- 
pear nearly round on the sky and never as flattened and fast-rotating 
systems. They cannot all be flat systems seen nearly face-on, as the 
observed fraction is too high (e.g. Tremblav & Merritt 1996). 



3.3 Influence of a dark-matter halo 

In the previous Section we derived the mass density for the dynam- 
ical models by deprojecting the observed stellar photometry under 
the simplifying assumption of a constant M/L. In the case of spiral 
galaxies, for which tracers of the potential are easy to observe, a 
number of observational and theoretical lines of evidence suggest 
that, either they a lso contain a substantial dark matter component 
(see lBinnevll2004 for a recent revie w), or that a modification to 
the N ewtonian gravity law is needed (Milgrom 19831: iBekensteinl 
l2004h . The situation is less clear regarding early-type galaxies, 
where it appears that in at least some galaxies dark matter is im - 
portant JCarollo et al]|l995t iRix et aljl?997t IGerhard et al1l200ll) 
while in others it may not be present a t all iRomanowsk v et alj 
l2003t iFerreras. Saha. & Williams] Eo05). In the simple case in 
which the dark matter density distribution follows the stellar den- 
sity, the adopted dynamical models will be correct, but the mea- 
sured M/L will include the contribution of the stellar component 
as well as the one of the dark matter. However in general the dark 
matter profile may not follow the stellar density, so the assumptions 
of our models could be fundamentally incorrect. 

We briefly investigate the effect of a possible dark matter halo 
on our results using a simple galaxy model. We assume the galaxy 
to b e spherical l y sym metric and the stellar density to be described 
by a lHemquistl 11990) profile (of total mass M = 1), which reason- 
ably well approximates the density of real early-type galaxies: 



p*(r) 



(3) 



27r r(r + l) :i ' 

Here the density goes as p* oc r _1 for r <C 1 and as p* oc r~ 4 for 
r ;§> 1. The dark matter contribution is represented by a logarithmic 
potential, 



ln(r 2 +ro) 



(4) 



which is the simplest potential producing a flat circular velocity vo 
at large radii (r ^> rn), as observed in real galaxies. 

With these assumptions, the projected second moment fi2 of an 
isotropic model, which is equal to the projected velocity dispersion 
a v due to th e spherical symmetry, is given, e.g., by equation (29) 
of lTremaine et alJl 19941) . Substituting the adopted expressions, we 
obtain: 



4(R) = 



ttTI{R) 
1 



+ 



4 



r(r + l) 2 (r 2 +r 2 ) 



\fi' 



R 2 



(r + 1) 3 



dr. 



(5) 



where I(R) is the stellar surface brightness at the projected radius 
R, T is the stellar M/L, and we assumed a gravitational constant 
G = 1. Thi s integral can be evaluated analytically (I(R) is given 
in Hernauist 199C), but it is trivially computed numerically. 

The results from this formula, evaluated for four different val- 
ues of Vo = 0, . . . , i>o,max, are shown in the bottom panel of Fig.|5| 
The characteristic radius ro = 3.5 of the dark halo and the value of 
vo,max = 0.45 were chosen so as to make the circular velocity v c 
of the stellar plus dark matter component as constant as possible in 
the range R — 1.815 — 5.446 (top panel of Fig.[5J, which corre- 
sponds to 1-3 half-light radii R c . This simple choice tries to mimic 
th e typical ob servations of real early-type galaxies (e.g. Figure 2 
of IGerhard et alj|200lh . Under these assumptions the dark matter 
represents 16% of the total mass inside R c and 52% inside 3R C . 
These values are ro ughly consistent with the values reported by 
Gerhard et al. (2001) for real galaxies. 
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Figure 5. Top Panel: black lines, from bottom to top, circular velocity v c , 
as a function of radius _R, of a spherical Hernquist galaxy model, with a 
dark halo described by a logarithmic potential, for four different values of 
the asymptotic halo velocity vq = 0, 0.15, 0.30, 0.45. The contribution to 
the circular velocity of the dark halo component alone is shown with the 
red lines. The three vertical dashed lines indicate the position of 1, 2 and 
3 half-light radii of the model. Middle Panel: Projected velocity dispersion 
o-p for the same four models shown in the top panel. The lowest curve is the 
Hernquist model without a dark halo (vo = 0). Bottom Panel: same as in 
the middle panel, for the luminosity-weighted aperture velocity dispersion 
<tr. This is the same quantity which is plotted in Fig.l2lfor the real galaxies. 



For this simple model, and assuming the largest dark matter 
contribution, the projected velocity dispersion a p is increased by 
22% at R a (the typical maximum radius sampled by our SAURON 
kinematics) and by 9% at Re/2, with respect to a model without 
dark halo. The luminosity-weighted velocity dispersion inside a cir- 
cular aperture or is increased by 8% inside R c and by 4% inside 
R c /2 (bottom panel of Fig.[3}. 

To get an indication of the overestimation of the stellar M/L 
we would obtain by fitting constant M/L models to galaxies with 
dark halos, we logarithmically rebinned the a p profile of the model 
with dark halo (to mimic the spatial sampling of the real data), and 
fitted to it the model without dark halo, in the radial range 0~R C . 



We obtained that the true stellar M / L would be overestimated by 
~ 12%, in the case of maximum dark halo. This shows that, under 
these hypotheses, the dark matter would start producing small but 
measurable variations in the global M / L. 



3.4 Three-integral dynamical modeling 

As we mentioned in the introduction, part of the observed FP tilt 
could come from non-homology in the orbital distribution of galax- 
ies, as a function of luminosity. To test if this is the case, we 
compared the M/L derived using the quicker and simpler two- 
integral leans models, with general stationary and axisymmetric 
three-integral dynamical models. 

The three-integral dynamical modeling we adopted is based 
on Schwarzschild's (1979) numerical orbit-superposition method, 
which is able to fit all kinematic and photometric observations 
jRichstone & Tremainelll988t iRix et aljfl997l: Ivan der Marel et alj 
1998). A similar approach was adopted by us and by other groups 



analyse the stellar orbital distribution (e. 


e. ICretton et al 


19991: 


CaoDellari et al. 20021: IVerolme et alJl2002i 


Gebhardt et al 


2003; 


VallurietalJl2004l 


iThomas et alJ [2004 


). The method provides a 



general description of axisymmetric galaxies, its stronger assump- 
tion being that of the stationarity of the galaxy potential. In practice 
the method only requires the potential not to vary on the time-scale 
required to sample the density distribution of an orbit. As the dy- 
namical time-scale in ellipticals is of the order of a few 10 7 yr at 
R c , while galaxy evolution generally happens on much longer time- 
scales, this assumption is expected to be valid in the regions that we 
sample with our kinematics. 

In our implementation of the Schwarzschild method we use 
the MGE parameterization to describe the stellar de nsity and to 
cal culate the gravitationa l potential, as in [Cappellari et al. (2002) 
an d lVerolmeet all 12002). However, we use here a new code which 
was specifically designed with the improved quality of our input 
integral-field kinematical data in mind. We verified that the new 
code produces the same results as the old one (within the nu- 
merical approximations of the method), when identical conditions 
are adopted. The key updates to the modeling method, with re- 
spect to the description given in Ivan der Marel et alj ll998t) and in 
ICappellarietai]<2002l) . consist of: 

(i) The orbit library required by the method has to be sampled 
by choosing the starting points of the orbits so as to cover a three 
dimensional space defined by the isolating integrals of motion, the 
energy E, the z-component of the angular momentum L z and the 
third integral 73. In previous works, at each E, we sampled orbits 
linearly in L z , and for each (E,L Z ) pair, an angle w as used to pa- 
rameterise I s (see Fig. 5 of lvan der Mare l et al. 1998). Here at each 
E we start orbits from the meridional plane with vr = v z = and 
L z > with positive initial coordinates arranged in a polar grid 
(linear in eccentric anomaly and in radius) going from R = z = 0, 
to the curve defined by the locus of the thin tube orbits (Fig. [6), 
to reduce duplications. The reason for our change is that we want 
to sample observable space as uniformly as possible, as opposed to 
sample the space of the integrals of motion, which is not directly ob- 
servable. With the new scheme the cusp singularities, which char- 
acterize the observables of orbits (see Fig. 1.2 of ICappellari et all 
2004) are distributed evenly on the sky plane. 

(ii) Ideally the orbit library in Schwarzschild's method should 
provide a complete set of basis functions for the DF. Individual or- 
bits corresponds to DFs that are 5-functions and an infinite num- 
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ber of them can indeed represent any smooth DF. However only a 
limited number of orbits can be used in real models, which then 
have a non-smooth DF and generate unrealistically non-smooth ob- 
servables. To alleviate this problem, instead of using single orbits 
as basis functions for the DF in the models, we use 'dithered' or- 
bital components, constructed by coadding a large bundle of orbits 
started from different but adjacent initial conditions. The starting 
points are arranged in such a way that the (regular) orbits differ 
in all three integrals of motion (Fig. [5}. This technique produces 
smoother components, which are also more representative of the 
space of integrals of motions they are meant to describe (a similar 
approach was used fo r spherical models bv lRichstone & Tremaind 
ll988tlRix et alll997l) . With dithered orbits we can construct mod- 
els which better approximate real galaxies than would otherwise be 
possible with single orbits, due to computational limitation. Dither- 
ing is crucial for this work, as otherwise the numerical noise gener- 
ated by the relatively small orbit number would limit the accuracy 
in the determination of the M/L. 

(iii) We do not use Fourier convolutions to model PSF effects. 
Given that the computation of the orbital observables is intrinsi- 
cally a Monte Carlo integration, we also deal with the PSF in the 
same way. The phase space coordinates are computed at equal time 
steps during the orbit integratio n (using the DOP853 routine by 
lHairer. Norsett. & Wanneilll993h . They are projected on the sky 
plane, generating the triplets of values (x',y',v' z ), which are re- 
garded as 'photons' coming from the galaxy. The sky coordinates 
(x' , y') of each photon are randomly perturbed, with probability de- 
scribed by the MGE PSF, before being recorded directly into any of 
the observational apertures. In this way no interpolation or interme- 
diate grids are required and the accuracy of the computation of the 
observables is only determined by the (large) number of photons 
generated. Voronoi 2D-binning is treated in the model in the same 
way as in the observed data (Section f2.3> . 

For the sampling of the orbit library in the present paper we 
adopted a grid of 21 energies, and at every energy we used 8 angu- 
lar and 7 radial sectors, with Aether = 6 (as in Fig.|6j. This means 
that the final galaxy model is made of 21 x (8 — 1) x 7 x 2 x 6 3 = 
444, 528 different orbits, which are bundled in groups of 6 3 = 216 
before the linear orbital superposition. From each orbital compo- 
nent we generate ~ 2 x 10 s 'photons' on the sky plane, which are 
used for an accurate Monte Carlo calculation of the observables in- 
side the Voronoi bins. At each bin position we fit V, a and four 
Gauss-Hermite parameters up to ha. 

Some of the models whose M/L we use in this paper were 
already presented elsewhere. An application of this new code to 
the modeling of the 2D-binned kinematics of the elliptical galaxy 
NGC 4473, which contains two counterrotating spheroids, and a 
study of the orbital distribution of the giant elliptical galaxy M87 
was presented in;Cap pellari & McDermidl2005l) . The modeling of 
the elliptical galaxy NGC 2974, and detailed tests of the ability of 
this orbital sampling sche me to recover a g iven M/L and a realistic 
DF are presented in Krainov Fc et aljl2005t) . 

In contrast to studies of spiral galaxies, where the inclination 
can be inferred from the observed geometry of the disk, there is 
usually no obvious way to infer the inclination of early-type galax- 
ies, and it therefore is one of the free parameters of the models. 
iKrainovic et alj 120051) showed that there is evidence for the incli- 
nation to be possibly degenerate using axisymmetric three-integral 
models, even if the LOSVD is known at every position on the pro- 
jected image of the galaxy on the sky. Here we confirmed this result 
using a larger sample of galaxies, in the sense that the models are 
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Figure 6. Visualization of our orbital starting scheme and dithering at a 
given energy E in the meridional plane (R, z) of an axisymmetric galaxy 
(cf. Fig. 5 of van der Marel et al. 1998). The large magenta dot marks the 
location of the circular orbit. At each E we start orbits from the meridional 
plane with vr = v z = and L z > 0, with initial coordinates arranged in 
a polar grid, linear in the eccentric anomaly and in radius (small red dots), 
going from R = z = 0, to the curve defined by the locus (green line) of the 
thin tube orbits, to reduce duplications. We then construct 'dithered' orbital 
bundles by grouping together (blue lines) every adjacent set of N^ itb or- 
bits and coadding them before computing the observables (at one E one can 
only visualize points or 6x6 in this plot). The first and last angular 

sectors have half size, to provide a uniform angular sampling across the R 
and z axes (considering the bi-symmetry of the model). Also shown is a rep- 
resentative regular non-resonant orbit, started near the centre of one orbital 
bundle, enclosed in its zero-velocity curve. The outermost black solid line 
is the equipotential, which includes every orbit at the given E. 



generally able to fit the observations within the measurement errors, 
for wide ranges of variation in inclination. For these reasons we fit- 
ted the M /L by constructing sequences of models with varying 
M/L, while keeping the inclination fixed to the best-fitting value 
provided by the leans models (Table 0. We show in Appendix A 
that the inclination derived from the leans models may be more ac- 
curate than just assuming all galaxies are edge-on, although we are 
aware it may not always be correct. In any case, we have shown in 
Fig.|4]that the results of this work do not depend critically on the 
assumed inclination. 

The Schwarzschild models use the same MGE parameters as 
the Jeans models (Section 13.21 and also include a central BH as 
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Figure 7. Comparison between the 25 M/ L determinations using the two- 
integral Jeans dynamical models and the three-integral Schwarzschild mod- 
els. An error of 6% in the model accuracy is required to explain the observed 
scatter along the best-fitting relation (solid line). The dotted line represents 
a one-to-one relation. 



predicted by the Mbh-o~ relation. More details on the modeling and 
the motivation of the changes, together with a full analysis of the or- 
bital distribution inferred from the models are given in a follow-up 
paper. The fact that the models can generally fit the observed den- 
sity and the SAURON integral-field kinematics of our galaxies with 
a constant M / L seems to confirm the suggestion of Section l3~3l ei- 
ther mass follows light in galaxy centres, or we cannot place useful 
constraints on the dark-matter halo from these data. This may be 
explained by the existence of an intrinsic degeneracy in the simul- 
taneous recovery of the potential and the orbital distribution from 
the ob served LOSVD, as discussed e.g. in Section 3 o flValluri et all 
12004), combined with the relatively limited spatial coverage. In 
a few galaxies the and ha Gauss-Hermite moments cannot be 
reproduced in detail by our models and may be improved e.g. by 
allowing the M/L to vary or extending to triaxial geometry. 



4 RESULTS 

4.1 Comparing M/L of two and three-integral models 

In this Section we compare the measurements of M / L obtained 
with the two-integral Jeans models described in Section l3~2l with 
the results obtained with the more general Schwarzschild models of 
Section l3~4l The usefulness of this comparison comes from the fact 
that the Jeans models are not general, but can be computed to nearly 
machine numerical accuracy, so they will generally provide biased 
results, but with negligible numerical noise. The Schwarzschild 
models are general, within the assumptions of stationarity and ax- 
isymmetry, but are affected by numerical noise and can be influ- 
enced by a number of implementation details. A comparison be- 
tween these two substantially different modeling methods will pro- 
vide a robust estimate of the uncertainties in the derived M/L, 
including systematic effects. The best-fitting correlation (Fig. 
shows a small systematic trend (M/L) Sc hw oc (Af/L)J^J 0,05 , 
while the scatter in the relation indicates an intrinsic error of 
6% in the models. Formal errors in the Schwarzschild M/L are 
generally not too dif ferent from this more robust estimate (e.g. 
iKrainovic et all2005l) . 

One can see in Fig.0that the galaxies with the highest M/L 



tend to show an (Af/L)s c hw which is systematically higher than 
(M/L) jo a n S (see also Table0. The difference in the M/L is likely 
due mainly to the fact that the Schwarzschild models use the full 
information on the LOSVD, while the Jeans models are restricted 
to the first two moments. For this reason, in the following Sections, 
unless otherwise indicated, we will use the term "dynamical M / L", 
or simply M/L, to refer to the (M/X)s c h w , which is expected to 
provide a less biased estimate of the true M/L. 

4.2 Observed correlations 

The data consists of N = 25 estimates of the dynamical M/L 
in the J-band and corresponding measurements of the luminosity- 
weighted second moment of the LOSVD <r c within the half-light 
radius. We fitted a linear relation of the form y = a + fix, with 
x — log(er/<7o) and y = log(Af/L), and minimizing 



E 

3=1 



(6) 



with symmetric errors Axj and A?/,, using the FITEXY routine 
taken from the IDL Astro-Library I Landsman 1993), which is based 
on a simila r routine bv |Press etall j 19921) . This is the same approach 
adopted by lTremdnee^lTl2002l) to measure the slope of the cor- 
relation between the mass Mbh of the central BH and <r c , and we 
refer to that paper for a comparison of the relative merits of alterna- 
tive methods. The M/L measurements are in versely proportional 
to the assumed distance, which is taken from iTon r^^^J (2001) 
and was adjusted to the new Cepheids zeropointof lFreedman^Hil] 
1200 ll) by subtracting 0.06 mag from the SBF distance moduli (Ta- 
ble Q}. For this reason the errors in the distance represent a lower 
limit to the errors in M/L. In the fit, we assigned to the M/L the 
distance errors quadratically coadded to the 6% modeling error de- 
rived in Section |4~T1 The a c values are measured from spectra of 
extremely high S/N, and have negligible statistical errors, but their 
accuracy is limited by template mismatch and possible calibration 
accuracy. We adopted an error in a e of 5%. An upper limit to the er- 
ror in a c can be derived by comparing our values with the central a c 
in published data (HyperLeda: 2 |Prugniel & S imien 1996), which is 
closely related to our cr c , but it is measured in a very different way. 
We found that errors of 7% in a are needed to explain the scatter 
with respect to a linear relation between log a c — log <r c . This scat- 
ter also includes very significant systematics due to the differences 
in the data, analysis and measurement method. The adopted 5% er- 
ror provides a conservative estimate of the true errors. Note that for 
th e gal axy M87 our <r c value is 20% smaller than the value adopted 
bv lTremaine et alJl2002l) . Our SAURON kinematics however agrees 
well with both the G -band and the Ca triplet stellar kinematics by 
van der Marel 1 1994). Adopting our smaller <r c M87 would become 
an outlier of the Mbh — c e correlation. 

We find a tight correlation (Fig.[8j with an observed rms scatter 
in M/L of 18%. The best-fitting relation is 



(M/L) = (3.80 ±0.14) 



(7) 



v200 kms _1 

The value of aq = 200 kms -1 was chosen, following 
iTremaine et alj J2002), to minimize the uncertainty in a and the 
correlation between a and (3. The best-fitting values and the un- 
certainties, here as in all the following correlations, were deter- 
mined after increasing the errors Ayj in the M/L by quadratically 

2 http://leda.univ-lyonl.fr/ 
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3 (l + /U,eV6), 



(8) 



log a c (km s" ) 



Figure 8. (M / L)-cr e correlation. The errors in the M/L values shown here 
and in the following plots represent the uncertainty in the distance as given 
bv lTonrvet ai]<200ll) . quadratically coadded to the 6% modeling error de- 
termined in Section l4*71 The errors on <r c are assumed to be 5%. The solid 
line is the correlation obtained by fitting all galaxies, while the dotted line is 
the fit obtained by excluding the five galaxies with log <r c < 2. The galaxy 
with the smallest <r c is M32. 



coadding a constant 'intrinsic' scatter so that (Ayj) 2 is replaced by 
(Aj/ i ) 2 + (Ay ) 2 ,tomakex 2 /^ = 1, where;/ = N-2 = 23 is the 
numb er of degrees of freedom in the linear fit fsee lTremaine et all 
2002, for a discussion of the approach). The constant intrinsic scat- 
ter implied by the observed scatter in this correlation is ~ 13% 
(Ay = 0.053). 

To test the robustness of the slope we also performed a fit 
only to the galaxies with log a c > 2, which is the range most 
densely sampled by FP studies (e.g. J0rgensen et al. 1996). The 
relation becomes in this case (M/L) oc a-*' 06 * 0,18 , which is 
steeper, but is still consistent with equation Q within the much 
larger uncertainty. The steepening of the correlation at high a c 
may also be due to a non perfect lin earity of the relation (e.g. 
Zaritskv, Gonzalez, & Zabludoff 2005), but this cannot be tested 
with our limited number of galaxies. The correlation computed us- 
ing the Jeans determinations is even tighter than equation Q, and 



gives (M/L) Jeans OC (7 C 



with an observed scatter of only 



15%. However, due to the possible bias of the Jeans M/L values, 
we will use the Schwarzschild determinations in what follows. 

The use of er determined from integral-field data and using a 
large aperture has the significant advantage that it is conceptually 
rigorous and accurate. However the (M/L)— a correlation does not 
appear to depend very strongly on the details by which a is deter- 
mined. We performed the same fit of Fig. [8] using as value of the 
dispersion the cr e /g measured from the SAURON data in an aper- 
ture of radius Re/8. In this case the best-fitting correlation has the 
form (M /L) oc cr°^ 5±0 ' 06 and the observed rms scatter in M/L is 
21%. As an extreme test we used as dispersion the inhomogeneous 
set of central values <r c , obtained from long-slit spectroscopic ob- 
servations, as given in the HyperLeda catalogue. In this case the 
best-fitting correlation becomes (M/L) oc o -0- 87 ± 07 an( j nas an 
observed scatter of 20%. Both correlations are consistent within the 
errors with the best fitting (M/L)-a e correlation, although they 
have a larger scatter. 

We also determined the correlation between the M/L and the 
secon d moment of the velo city, corrected with equation (18) of 
Ivan der Marel & Franxl 1 19931) 



to include the contribution to the second moment of a nonzero hi, c 
parameter. Equation {5} was obtained by integrating the LOSVD 
to infinite velocities. In practice the second moment of a param- 
eterised LOSVD with dispersion a and h\ = 0.1 is equal to 
a — (1, 1.11, 1.21) x a if one sets the LOSVD to zero for ve- 
locities larger than | V| > (2.45, 3,4) x a respectively. This shows 
that the correction has to be used with care, being highly sensitive 
to the details of the LOSVD at large velocities, where the LOSVD 
cannot be determined reliably. In fact a measured positive /14 pa- 
rameter, which means the LOSVD is narrower than a Gaussian at 
small velocities, does not necessarily imply that the wings of the 
LOSVD are precisely as described by the Gauss-Hermite parame- 
terization. We find a correlation consistent, within the errors, with 
equation Q and with comparable scatter. 

We also studied the correlation of the dynamical M/L with 
the galaxy luminosity (Fig. [5}- For this we used our own total mag- 
nitudes in the /-band It given in Table Q transformed into abso- 
lute magnitudes using the SBF distances and the galactic absorp- 
tion iSchl egel etaiT l998) given in Table Q The It measurements 
are dominated by possible systematic errors due to the extrapolation 
involved and the possible contamination by nearby objects. To ob- 
tain a realistic estimate o f the errors we com pared our It with the 
same value 3 measured by T onrv et afl feoOl) in the same photomet- 
ric band. A fit to the two sets of values showed no systematic trend 
with magnitude and a zeropoint shift of 0.09 mag. The observed 
scatter is consistent with errors of 13% in each measurement, and 
we assigned this error to our It values. A stringent internal test 
on the errors in the It measurements comes from the comparisons 
with the total luminosities obtained by coadding all the Gaussians 
of the MGE models of Section l3~T1 A correlation of the two quanti- 
ties shows no deviant points and no systematic trend in either slope 
or zeropoint, and is fully explained by an 8% random error in both 
quantities. This estimate is simi lar to the typ ical measurement er- 
rors of 9% reported bv lj0rgensen et al] ll996l) . We also compared 
It with the corresponding v alues of Bt in the _B-ba nd as given 
by the HyperLeda catalogue iPrugniel & Simien 1996). We found 
a smooth variation of the B — I colour with luminosity, consistent 
with the quoted measurements errors in this quantity and indicating 
the absence of significant systematic errors or any deviant point. 
The fit was performed in the same way as above and provided the 
relation: 



(M/L) = (2.35 ±0.19) 



Li 



10 10 ^, Q 



(9) 



This correlation has an observed rms scatter in M/L of 3 1 % which 
is much larger than that of the correlation 0. In particular there is 
a significant scatter for the eight galaxies in the luminosity range 
—21.5 < Mi < — 19.5 mag. Given the fact that all galaxies in the 
SAURON sample were selected to have Mb < — 18 mag, which 
translates in the J-band into Mi < — 20 mag, the fit depends very 
sensitively on the single M/L of M32, which also has a small dis- 
tance error. Moreover there is some indication that the dwarf E 
galaxies like M32 may popula te a different FP from t he normal E 
galaxies iBurstein et alJI 19971 : but see iGraham & Guzmarj fcoOS). 
For these reason we excluded M32 from this fit and all the follow- 
ing fits involving luminosity or mass. 



It = mi 



Ni in the notations of Table 1 of lTonrv et al]<200lh 
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Figure 9. (M/ L)-Lj correlation. Note that all galaxies, with the exception 
of M32 are brighter than the limit Mj ss — 20 mag, in absolute magnitudes, 
which is set by the luminosity selection criterion of the SAURON survey. Due 
to the sensitivity of the fit to the single M/L of M32, this galaxy was not 
included in this and in the following fits involving luminosity or mass, but it 
is included in the plots, to allow a comparison with Fig. 151 



4.3 Understanding the tightness of the correlations 

The decrease of the scatter in M/L, when going from the (M/L)- 
L correlation to the (M/L)-a one is reminiscent of what happens 
to the similar correlations involving the BH mass Mbh instead of 
the M/L. This decrease of the scatter shows that at a given galaxy 
luminosity the M/L is higher when the a is higher. One may ap- 
ply the argument s used to understand the tightn ess of the Mbh- 
a relation iGebhardt et al. 2000: iFerrarese & Merritd l2000) to the 
(M/L)-a relation as well. There can be two main reasons for the 
(M/L)-a relation to be tighter than the (M / L)-L relation: (i) the 
M/L may correlate with galaxy mass, but the luminosity may not 
be a good mass estimator, due to differences of M/L at a given lu- 
minosity; in this case variations of a would reflect variations of 
M/L; (ii) the M/L may be related to the galaxy compactness, 
which implies a higher a at a given mass. I n the case of the Mbh- 
a rela tion both effects seem to play a role: Kormendv & Gebhardt 
1200 ll) showed that the gala xies with higher a are indeed mo re com- 
pact than average (see also Graham, Truiillo, & Caon 2001). Subse- 
1 120031) 



quently Marconi & Hun] 120031) and lHaring & Rixl 120041) showed 
that much tighter correlations for Mbh are obtained using the 
galaxy mass instead of the luminosity. 

It is clear from Fig.|5|that there is a spread in M/L at a given 
galaxy luminosity, so the first argument above must play a role in 
tightening the correlation. To test the importance of this effect we 
verify how much better the M/L correlates with galaxy mass after 
multiplying each luminosity by its measured M / L. The best-fitting 
correlation (Fig. llQI is given by 



(M/L) = (1.78 ±0.16) 



M 



10 10 A% 



(10) 



This correlation has an observed scatter in M/L of 24%. It is tighter 
than the (M / ' L)-L correlation, but still significantly less tight than 
the (M/L)-a one. 

An alternative way to verify the importance of the M/L vari- 
ations in increasing the scatter of the (M/L)~L relation is to use 
the galaxy luminosity in the near infrared, instead of in the /-band. 
The near infrared better samples the light emitted by the bulk of the 
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Figure 10. (M/L)-M correlation. Here the /-band luminosity of each 
galaxy was converted into mass by multiplication by its M/L. The least 
massive galaxy is M32, and was not included in the fit. 



0.8 



I 0.6 



0.4 



c 



0.2 



0.0 




-26 



Figure 11. (M/ L)-Lk correlation. Same as in Fig.|9| but here the M/L is 
plotted against the 2MASS galaxy absolute magnitude in the X-band. The 
faintest galaxy is M32 which was not included in the fit. 



galaxies' old stellar population. Moreover in that wavelength range 
the M /L variations with mass are smaller (e.g. Bell & de Jong 
|200jh. These are apparently the reasons why iMarconi & Hunt 
12003) found a much tighter Mbh-L correlation in the /C-band 
than in the _B-band. 

We want to see if the same tightening happens in our case with 
the M/L. The correlation between the /-band M/L and the K- 
band luminosity, extracted from the 2MASS Extended Source Cat- 
alog (XSC; Jarre tt et alj2000l) . is shown in Fig.[U] We adopted the 
same constant error of 13% in the luminosity, as we used in the 
/-band, instead of the errors given in the XSC, which have an un- 
realistically small median value of 2%. The best fitting relation is 
nearly identical to equation j5J having the form 

0.32±0.O5 

(M/L) = (1.88 ± 0.20) f T7 ^ ) , (ID 



adopting M K q = 3.28 mag. This relation has indeed a slightly 
smaller rms scatter in M/ L, of 28%, than the /-band correlation. 

Finally, to test the importance of the second argument above, 
namely the effect of the galaxy compactness in tightening our 
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Figure 12. (M/ L)-Lj correlation as in Fig. |5] The radius of the circles, 
centered around each measurement, is proportional to the galaxy R c in pc. 



(M/L)~a correlation, we show in Fig. ll2l the same (M/L)-L re- 
lation as in Fig. 00 but we also overplot the effective radii R c of the 
galaxies (in pc). The eight galaxies with — 21.5 < Mj < — 19.5 dif- 
fer substantially in their R c . Indeed the galaxies with higher M/L 
are also the most compact at a gi ven luminosity. This explains w hy 
they also have the highest a (see lKormendv & Ge bhardt 20(3, f° r 
a discussion), and it shows why, at least for our limited galaxy sam- 
ple, the (M/L)-a correlation is even tighter than the (M/L)-M 
relation. 

To summarize the results of this Section, we found that (i) the 
tightest correlation is the (M/L)-a relation, followed by (A'l/L)- 
M and (M /L)-L K - The (M/L)-Li relation is the least tight. This 
sequence of increasing scatter suggests that the M / L correlates pri- 
marily with galaxy mass and, at given mass, with galaxy compact- 
ness. However we caution the reader that the size of our sample is 
too small to draw definitive conclusions. In fact, with the exception 
of the (M/L)-a relation, the difference in scatter between the dif- 
ferent correlations is only at a few sigma level. A larger sample is 
needed to make this result stronger. 



4.4 Understanding the tilt of the FP 



The existence of the FP of the form R c oc a a Ig, combined with 
the homology and virial equilibrium assumption M oc a 2 R c and 
the geometric definit ion L oc I e R 2 , yields a pr ediction for the 
(M/L) F p iDressler et alJl987tlFaber et alll987l) . Here we adopt 
a FP of the form: 



Re 



(12) 



as determined in the Gunn r-band by Jorge nsen et all II 996) from a 
sample of 225 early-type galaxies in nearby clusters. We choose 
these numerical values because the photometric parameters R c 
and I a for this FP were measured from a homogeneous sample 
of E and SO galaxies, in the same way as we did (Section 12.21 . 
This FP determination agrees within the errors with other deter- 
minations in similar photometric bands (see Colless et al. 2001, 
and refere nces therein) and also with the /-band determination by 
IScodeggio. Giovanelli. & HavneJ ll997t) . so the choice of the FP 
coefficients is not critical. One impo rtant exception is the FP de- 
termination bv lBernardi et alj 120031) . derived from a large sample 
of galaxies from the Sloan Digital Sky Survey. The parameters of 



that FP are inconsistent with most previous determinations, and the 
reason for this discrepancy is currently unclear. Considering also 
that their photometric parameters were derived by fitting the galaxy 
profiles and not from a growth curve, as in this and other FP works, 
we decid ed not to consider their values here. Taking the numerical 
values bv ljetrgensen et all il996). one can derive the following re- 
lations for the (M/L)pp, as a function of the combination R e a or 
equivalently galaxy mass, and galaxy luminosity : 



(M/L) 



FP OC 



R, 



-1-1/0 a 2 + a/f3 



D 0.22±0.03 0.49±0.09 

R n a 



(M/L)j. 



_ R -(2+ Q )/(2/3)-2 M l+ Q /(2/3) 
^,-0.O2±O.06jy^-0.24±0.05 
j-l/2-(14-2,3)/c«£-l/2+l/<* 
j0.02±0.04^0.31±0.05 



(13) 
(14) 



The scatter in the relations amounts to 23% between M/L and 
R°- 22 a 0Ag or mass, and 31% between M/L and L. Within the er- 
rors the (A//L)fp can be written as a simple power-law of either 
luminosity or mass alone. This simple dependence is what is usu- 
ally known as the 'tilt' of the FP. The (A//L)fp cannot be written 
explicitly as a simple function of cr, however [j0rgensen et aljfT996l) 
note that "the inclusion of R c is not even essential" as a direct cor- 
relation between (M/L)pp and a gives a correlation: 



(M/L)p P oc cr° 



(15) 



with an observed scatter virtually unchanged of 25%. 

As mentioned before, the correctness of the FP prediction for 
the M/L depends on the homology assumption, regarding both the 
surface brightness and the orbital distribution of early-type galaxies 
(e.g. lCiotti et alj|l996t) . Our dynamical M/L is based on detailed 
photometric models, so it is independent of the spatial homology 
assumption. Moreover we include the full LOSVD to estimate the 
M/L and we do not need to make any extra correction for galaxy 
rotation. Finally the good agreement between the M/L derived 
with two-integral Jeans and three-integral Schwarzschild models 
shows that orbital non-homology does not have an important effect 
on the M/L fSection f4. H . We are thus in the position to compare 
directly our M / L determinations, which do not depend on any ho- 
mology assumption, with the homology-dependent (M/L)fp. 

Comparison between the FP tilt given by equation 1 141 and 
relation j5J shows that the dependence of the dynamical M/L on 
luminosity is the same. Taking the measurement errors into account 
this indicates that non-homology cannot produce more than ~ 15% 
of the observed tilt. For the more reliable correlations involving a, 
comparison of the virial prediction of equation I15i and the model- 
ing result of equation Q indicates that non-homology can account 
for at most ~ 7% of the tilt. Both relations consistently imply that 
the FP tilt reflects essentially a real variation of the total M/L in 
the central regions of galaxies, which can be due to variations in 
the galaxies' stellar population and/or in the dark matter fraction. A 
similar conclusion was reached bv lLanzoni et alJ l2004ll from gen- 
eral considerations about the observed scaling relations of early- 
type galaxies. The comparison of this Section however involves un- 
certainties due to the fact that the tilt may depend on the sample 
selection criteria. To assess if this plays an important effect we per- 
form in the next Section a direct comparison of the virial predictions 
for the M/L derived from our own galaxy sample. 
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can provide a direct 'calibration' of the virial mass, and M / L esti- 
mator (which are often used only in a relative sense): 



0.2 0.4 

log (M/L) vlr 

Figure 13. Comparison between the (M/L) v i T = /3iJ c cr^/(L G) derived 
from the virial assumption and the M/L obtained from the Schwarzschild 
models. The values of (M/L) v i v were scaled to match the dynamical M / L, 
and the best-fitting factor is f) = 5.0 ± 0.1. The solid line is a fit between 
the two quantities, while the dotted line represents a one-to-one correlation. 



4.5 Comparison with virial predictions of M/L 

An alternative way to test the validity of the virial and homology 
assumptions and their influence on the FP tilt is to compute the 
'observed' virial (M/L) v i T oc R e ag/L and to compare it directly 
to the M/L derived from the dynamical models. This has the ad- 
vantage that it can be performed on our own galaxies and does not 
involve any choice of FP parameters or selection effects. We fitted 
the correlations of (A'//L) v i r , in the /-band, with a and with lumi- 
nosity, obtaining: 



(M/L) viT oc a c ° 
(M/L) vir oc L° 



(16) 
(17) 



Equations 1 161 171 have an observed rms scatter of 19% and 27%, 
and are fully consistent with the FP determinations in equations 1 1 5i 
and 1 141 respectively. The scatter in these correlations of virial de- 
terminations is comparable to the scatter derived using the full dy- 
namical models and, as in that case, appears dominated by the in- 
trinsic scatter in M / L. 

Finally, the most direct way of measuring the accuracy of the 
homology assumption is to compare (M/L) v i r with the M/L from 
the dynamical models. The correlation is shown in Fig. ll3l and has 
the form 



(M/L) oc {M/L)lf T s±om . 



(18) 



The observed slope is consistent with the determination based on 
the FP (Section l44l . implying that both the structural and the orbital 
non-homology contribution cannot represent more than ~ 15% 
of the FP tilt (neglecting possible selection effects in our sample, 
which are very difficult to estimate). Contrary to all correlations 
shown before, the scatter in this correlation is not influenced by the 
errors in the distance, as both M / L estimates use the same distance. 
Adopting an intrinsic accuracy of 6% in the M/L determinations 
(Section f4.lt . the scatter in (Af/L) v i r required to make \ 2 / v = 1 
is 14%. The correlation between the virial and the Jeans estimates 
gives (M/L)j, 
the correlation 1 1 8i 

Comparing the virial and Schwarzschild M / L estimates we 



(M/L)vir = 13 R e a 2 J(LG). 



(19) 



oc (M/L)° v f 4±0M with very similar scatter to 



The best fitting scaling factor is f3 — 5.0±0.1. We can compare this 
value with the predictions from simple dynamical models, as done 
by a number of previous authors (e.g. Michard 1980). For this we 
computed the theoretical predictions for f3 from spherical isotropic 
models described by the Sersic profile R 1 ^", for different values 
of n. The computation was performed using high-accuracy MGE 
fits to the Sersic profiles, obtained with the routines of lCappellaril 
12002). From the fitted MGE models, which reproduce the pro- 
files to better than 0.05%, the projected a values can be computed 
with a single one-dimensional numerical integration. The projected 
luminosity-weighted a was then integrated within a circular aper- 
ture of radius R c to compute a c which is needed to determine the 
scaling factor [3. In the range n — 2-10 the predicted /3 parameter 
is approximated to better than 3% by the expression 

(3{n) = 8.87- 0.831 n + 0.0241 n. (20) 

(cf. iBertin. Ciotti. & Del Principell2002h . The precise value pre- 
dicted for a R 1 / 4 profile is j3 = 5.953 (the value becomes = 
5.872 with a BH of 0. 14% of the galaxy mass as in lHaring & Rixl 
2004), while the observed value of j3 ~ 5.0 would correspond to 
a Sersic index n ~ 5.5. However the predictions of equation 1201 
only apply in an idealized situation and do not take into account the 
details in which (M/L) v ir is measured from real data and the fact 
that galaxies are not simple one-component isotropic spherical sys- 
tems. From our extended photometry (Section l2~2l we measured the 
n values for the 25 galaxies of our sample by fitting the observed 
radial surface brightness profiles. The derived Sersic indices span 
the whole range n — 2-10 and will be presented in a future paper. 
From the observed variation in the profiles (5 should be expected to 
vary by a factor ~ 2.5 according to the idealized spherical model. 
In practice we find no significant correlation (linear correlation co- 
efficient r ~ —0.13) between the measured (3 (the value required 
to make (M/L)s c hw = PR c a^/(LG)) and the one predicted by 
equation <20> . This shows that the idealized model is not a useful 
representation of reality and cannot be used to try to improve the 
(M/L)vir estimates. An investigation of the interesting question of 
why the /3 parameter appears so constant in real galaxies and with 
realistic observing conditions goes beyond the scope of the present 
paper. 

The results of this Section show that the simple virial estimate 
of M/ L, and correspondingly of galaxy mass, is virtually unbiased, 
in the sense that it produces estimates that follow a nearly one-to- 
one correlation with the M/L computed from much more 'expen- 
sive' dynamical models. This result has implications for high red- 
shift studies, as it shows e.g. that one can confidently use the virial 
estimator to determine masses at high redshift, where it is gener- 
ally unfeasible to construct full dynamical models. It is important 
to emphasize however that this result strictly applies to virial mea- 
surements derived as we do in this paper, namely using 'classic' 
determination of R c and L via R 1/4 growth curves, and with <r e 
values measured in a large aperture. 



4.6 Comparison with previous results 

To understand whether the tight (M/L)~a c correlation is really due 
to the improved quality of the kinematical data and photometric 
data, we determined the observed scatter in the Johnson i?-band 



© 2006 RAS, MNRAS OOO.rflEsl 



The SAURON project -IV 15 




o.o r i i i , , i , , , i 

1.8 2.0 2.2 2.4 2.6 

log g c (km s" 1 ) 

Figure 14. Same as in Fig. |8] for the M/L measurements obtained by 
Ivan der Marel 1 1991) on a different galaxy sample. His R-band M/L val- 
ues have been rescaled to the same SBF distances and zeropoint we use here 
and converted to the /-band assuming {M / L) R / {M / L) I 1.27. The 
solid line is the best-fitting relation to the data, while the dotted line shows 
for comparison relation 0. 

M/L determinations by Ivan der Marell 1 199 ll) . who used a simi- 
lar Jeans modeling method, on a different sample of bright ellip- 
tical galaxies (five in common with our sample). From his sam- 
ple of_37_galaxies_we extracted the 26 for which SBF distances 
bv iTonrv et alj J2001) exist and we performed a fit as done in 
Fig. [8] We used in the correlations his mean velocity dispersion 
a and his "improved" (M/ L) % ^ p (his column [14] from Table 1 
and column [18] of Table 2 respectively), which we rescaled to 
the SBF distances, and converted to the /-band assuming an R- 
band absolute magnitude of the Sun of M R q = 4.42 mag and 
a typical colour for ellipticals of R — I ~ 0.6, which implies 
(M/L)r/(M/L)i w 1.27. As in Section l4~2l we assumed a 5% 
error in a, and assigned to the M / L the errors in the SBF distance, 
quadratically coadded to a constant error as to make x 2 /v — 1. The 
best-fitting relation is shown in Fig. ll4l and has the form: 

/ a \0.79±0.17 

(M/L) = (3.54 ± 0.17) ( * A (21) 
V 200 km s / 

which is in good agreement with our equation 0, but has much 
larger uncertainties and an increased observed scatter of 28%. This 
increased scatter must be attributed mainl y to th e differences in the 
quality of the data used by van der Marel 1 1991) and in the present 
paper. For the five galaxies (NGC 3379, NGC 3608, NGC 4374, 
NGC 4486, NGC 5846) in common with this paper we measured 
relative differences e = 1 - (M /L) vdM /(M /X)c ap of 26%, 14%, 
-11%, -1%, -27% respectively. 

4.7 Comparison with stellar population (M/L) pop 

We have established that the dynamical M/L can be accurately 
determined from our data and models, and that a tight correlation 
exists with <r c , which also explains the appearance of the FP. The 
next obvious question is whether the observed M/L variations are 
mainly due to a change in the stellar population or to differences in 
the dark matter fraction in galaxies. A number of authors have ad- 
dressed this question usually using simple dynamical models com- 
bined with stellar population models. Sometimes these investiga- 
tions have produced outcomes t hat are not entirely co nsistent be- 
tween each other ("e.g. lKauffmann et alJl2003t lPadmanabha n et alj 



2004: lDrorv. Bender. & HoprJ2004l) . Here we use accurate dynam- 
ical M/L and (M/L) pop estimates from singl e stellar population 
(SSP) models, similarly to Gerhard et al. 1 2 00 ll) . 

Using the (M/ L) pop predictions for the SSP models 4 by 
IVazdekis et al] j 19 9(1 here after VZ96), and adopting as initial ref- 
erence the lsalpeterl il955t) stellar initial mass function (IMF), we 
derived (M/L) Do p from the SAURON Ms 6, Fe5015 and H/3 line- 
strength indices iKuntschner et al. 2006, Paper VI). To minimize 
the uncertainties in the absolute calibration of the line-strength in- 
dices to model predictions, we determined the index predictions for 
different SSP ages and metallicities by measuring them directly on 
the flux calibrated model spectral energy distributions produced by 
VZ99 (see discussion in Kuntschner et al. 2002). We hereby cir- 
cumvent the use of the so-called Lick fitting functions which can 
only be used in conjunction with uncertain offsets accounting for 
differences in the flux calibration between models and observations. 
The model spectra of VZ99 were broadened to the Lick resolution 
using the same procedure we adopted for the SAURON spectra. In 
order to minimize the influence of non-solar abundance ratios we 
derived the (M / L) pop estimates from a [MgFe50] vs H/3 diagram 
(Fig-EJ. The index [MgFe50] = (0.45 x Mg6 + Fe5015)/2 
is defined such that it is largely insensitive to abundan ce ratio 
variations, while H/3 is only wea kly dependent on it fTrager et all 
l2000r,lThomas. Maraston. & Kornl2004 . For each galaxy, the line- 
strength measurements were extracted from the same high signal- 
to-noise SAURON spectrum from which the er c were derived, i.e. by 
luminosity-weighting all the spectra within R c (Section f2.3> . 

To test the robustness of the (M/L) pop predic- 
tions we repeated the ana ly sis o f Fig. 1151 using the mod- 
els of iBnizual & Charlol 120031 BC03), and those by 
j Thomas. Maraston. & Bended 120031 TMB03) and iMarastorJ 
12005), adopting the Salpeter IMF. The absolute values of the 
(M/L) pop agree to within 20% between the different models over 
the line-strength index range relevant for our galaxies. However 
the detailed agreement is better between the TMB03 and VZ96 
models, while the predictions of BC03 tend to be noisier. In this 
paper we choose to use the VZ96 and VZ99 models because for 
them predictions of the spectral energy distribution at moderate res- 
olution are also available which is helpful to minimize calibration 
issues. 

We find that for single-burst stellar populations, the contour 
levels of constant (M/L) pop in Fig. ll5l tend to be nearly horizontal, 
in the sense that (M/L) pop is essentially a function of H/3 alone. 
This implies that, under the model assumptions, if the variations in 
the dynamical M/L are driven by the variation in the stellar popu- 
lation, a good correlation should exist between M/L and H/3. This 
is tested in Fig. l 161 where we plot the measured M/L and H/3 mea- 
surements and compare them to the (M/L) pop versus H/3 predic- 
tions, for a set of models with a large spread in age (2-17 Gyr) and 
metallicity ([M/H] = [-0.68,-0.38,0.00,0.20]). As expected 
from Fig. ll5l the envelope of the different model curves traces a tight 
relation, nearly linear in logarithmic coordinates (this is less true if 
one also considers models with more extreme [M/H] values). The 

4 Available from http://www.iac.es/galeria/vazdekis/ We used the model 
predictions for the (M / L) pop , that were updated in February 2005 to take 
into account the contribution of the remnants and the mass loss during the 
latest phases of the stellar evolution. The adopted lower and upper mass 
cutoffs for the IMF are 0.01 and 120 Mq respectively, whereas the faintest 
star is assumed to be 0.09 Mq. 

5 The values were taken from Table 5 of Paper VI after removing the ap- 
plied Lick/IDS offsets listed in Table 3 of that paper. 
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Figure 15. Observed values of the line strength index log[MgFc50] ver- 
sus log(H/3) for our galaxy sample. Overplotted are the contours show- 
ing the corresponding predictions of log(Af/L) pop for the SSP models 
by VZ96 and VZ99, using the Salpeter IMF. The small crosses indicate 
the location where the model predictions are available, and from which 
the contours were linearly interpolated. Note the concentration of galax- 
ies near the contour level log(A//L) pop = 0.6. The contours appear 
nearly horizontal in this plot, and the same is true using the Kroupa 1 200 1 ) 
IMF. The two galaxies with the H/3 value lower than any model prediction 
(M87 and NGC 5846) may be contaminated by residual gas emission. Their 
(M / L) pop was estimated after increasing H/3 to the level of the model 
predictions. 

model relation follows the same trend as the data, suggesting that 
the variation of the stellar population is indeed an important fac- 
tor in determining the observed variations of the M/L. The span 
of M/L from 1-6 mainly reflects differences in the luminosity- 
weighted age of the stellar population. T he observed age tr end is 
consistent with evidences of downsizing iCowie et al. 1996) from 
high redshift studies (e.g.|di Serego Alighieri et al l2005l : lTreu et alJ 
2005: lvanderWeletall2005h . 

The accuracy of our M/L determinations allows us to go be- 
yond the general agreement, to detect significant deviations, and to 
exclude a simple one-to-one relation between the total M/L given 
by the dynamical models and the (M/L) pop predicted by the stel- 
lar population models. In particular, for some galaxies, the mea- 
sured total M/L is lower than any of the model predictions with 
the Salpeter IMF. For this reason, with the above caveats about the 
absolute scaling of the (M/L) pop values, and if we assume the 
IMF to be the same for all galaxies, we have to reject the Salpeter 
IMF (and any IMF with larger slope) as it gives unphysical results. 
lKrouDdl200ll) has measured the IMF and constrained the shape be- 
low one solar mass with the result that there are fewer low mass 
stars than indicated by the Salpeter law. The effect of using the 
Kroupa IMF is just to decrease all the (M / L) pop model predic- 
tions by ~ 30% (Alog(M/L) ~ 0.16). With this IMF, none of the 
dynamical M/L values lie below the relation defined by the model 
predictions, within the measurement errors. Although we cannot in- 
dependently confirm the correctness of this IMF we adopt it for the 
purposes of the following discussion. 

A more direct, but more model dependent, way of looking 
at the relation between the stellar population and the M/L vari- 
ations, is to compare M/L with (M/L) pop predicted using the 
Kroupa IMF (Fig. 1171 . The similarity of this figure with Fig. 1161 
which was obtained from the data alone, shows that the structure 
seen in the plot is robust, and does not come from subtle details 
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Figure 16. Observed values of the dynamical M/L versus the observed 
line-strength index log(H/3). There is a clear inverse correlation between 
the two quantities, however the scatter is significant, and cannot be ex- 
plained by measurement errors alone. The red squares and the blue dia- 
monds indicate slow and fast rotating galaxies respectively, as defined in 
Table IT1 The coloured lines show the predictions of log(Af/L) pop ver- 
sus log(H/3) from the SSP models of VZ96 and VZ99, using a Salpeter 
IMF. The blue, green and red solid lines correspond to a model age of 2, 7 
and 17 Gyr respectively. The dashed lines indicate a metallicity [M/H] = 
[-0.68,-0.38,0.00,0.20], the red dashed line being solar metallicity 
[M/H] = 0. The difference in the labelling of the left and right axes is 
meant to emphasize the fact that along the vertical axis we are comparing 
two different quantities: the dynamical or total M/ L, and the (M/ L) pop of 
the stellar population alone. All the population model lines essentially over- 
lap and define a tight relation between the two quantities, which is generally 
consistent with the measured trend. However the fact that some of the mea- 
sured total M/ L lie below all of the population model predictions shows 
that the adopted IMF is unphysical. Adopting a Kroupa IMF the (M/ L) pop 
model predictions decrease by ~ 30% (A log(Af/L) ~ 0.16) and all the 
observed M/L lie above or on top of the stellar population models. This 
indicates that a Kroupa IMF is consistent with the observations for all the 
galaxies in our sample. 

in the SSP predictions. Again the main result is a general correla- 
tion between t he total M/L and the stellar population (M / L) pop , 
consistent with lGerhard et all 1200 lL see their Fig. 14), albeit with 
a smaller scatter. The relatively small scatter in the correlation indi- 
cates that the IMF of the stellar population varies little among differ- 
ent galaxiej J _coiisistent with the results obtained for spiral galaxies 
bv lBell & de Jonsl <200ll) . All galaxies have (M / L) pop < (M/L) 
within the errors, but the total and stellar M/L clearly do not fol- 
low a one-to-one relation (green thick line). Unless this difference 
between M/L and (M/L) pop is related to IMF variations between 
galaxies, it can only be caused by a change of dark matter frac- 
tion within 6 one R c . The inferred median dark matter fraction is 
29% of the tot al mass, broadly consistent previou s findings from 
dynamics (e.g. Gerhard et al. 2001; Thomas et al. 2005) and grav- 
itational lensing (e.g. lTreu & Koopm anlll2l)o'4llRusiri*& Kochanek 
2005), implying that early-type galaxi es tend to have 'minimal ha- 
los' as spiral galaxies (Bell & de Jong 2001). The large number of 
galaxies in our sample with a similar value of log(A//L) pop ~ 0.5, 
and the clear evidence for a variation in the dark matte r fraction in 
those galaxies is consistent with the result of Padmanabhan et al. 

6 This statement is not entirely rigorous, in fact, as shown in Section l3*3l 
the presence of a dark matter halo can increase the M/L measured within 
R e in a way that is not simply related to the amount of dark matter within 
that radius. 
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Figure 17. Dynamical (i.e. total) M/L from the Schwarzschild modeling 
as a f unction of (M / L) POD using the SSP models of VZ96 and VZ99, 
with a iKroupJ 1200 1) IMF. The red squares and the blue diamonds in- 
dicate the slow and fast rotating galaxies respectively, as defined in Ta- 
ble □ The thick green line indicates the one-to-one relation. All galaxies 
have (M/L) pop < (M/L) within the errors, but the total and stellar M/L 
clearly do not follow a one-to-one relation. Dark matter is needed to explain 
the differences in M/L (if the IMF is not varying). The magenta arrows 
show the variation in the estimated (M / L) pop for the youngest galaxies 
(luminosity-weighted age < 7 Gyr), if a two-population model is assumed 
(see text for details). The (M / L) pop of the young galaxies would move 
closer to the one-to-one relation. Adopting the Salpeter IMF all the values 
of (M / L) p0 p would increase by Alog(M/L) ~ 0.16. This can be visu- 
alized by shifting the one-to-one relation to the position of the dashed line. 
In this case a number of galaxies would have (M / L) pop > (M/L) and 
this implies that the Salpeter IMF is unphysical (consistent with Fig. 1161 . 

120041) . However our results are not consistent with a constancy of 
the (M / L) pop for all galaxies, and a simple dark-matter variation, 
as we clearly detect a small number of young galaxies with low 
(M/L) pop and a correspondingly low dynamical M/L. Adopting 
the Salpeter IMF all the values of (M / L) pop would increase by 
Alog(A//L) ~ 0.16. This can be visualized by moving the one-to- 
one relation to the position of the dashed line. In this case a number 
of galaxies would have (M / L) pop > (M/L), implying that the 
Salpeter IMF is unphysical (as we inferred from Fig. 1161 . 

This analysis has been carried out in the context of SSP mod- 
els. We now consider a few caveats arising from this methodol- 
ogy. Low-level secondary star formation can alter the values of 
the (M/L) pop derived: essentially we underestimate (M/L) pop 
for these galaxies. This effect is expected to influence mainly the 
youngest galaxies, with the lowest (M / L) pop . To estimate the im- 
portance of this effect we experimented with combinations of two 
SSP models, considering a young component on top of an older 
solar-metallicity one. In Fig.[T71 

we indicate with an arrow the vari- 
ation in the estimated (M /L) pop , for the five youngest galaxies 
(luminosity weighted age less than 7 Gyr), assuming the stellar pop- 
ulation is composed by 90% (in mass) by an old 12 Gyr population 
with solar metallicity and by 10% by a younger population, whose 
age and metallicity is allowed to vary to reproduce the observed 
line strengths. We find that the artificial underestimation of the low- 
est (M/L) pop values can easily explain why these galaxies appear 
to have a high fraction of dark matter, contrary to the general trend. 
The adopted models use element abundances in the solar ratios, yet 
we know that many of our sample galaxies have enhanced ratios of 
Mg/Fe. Unfortunately, there are no models available which predict 



(M / L) pop as a function of abundance ratios. A possible influence 
on our results remains to be explored. 

Given our assumptions about the IMF, we conclude the fol- 
lowing: (i) There is a reasonable correlation between the dynamical 
and stellar population estimates of M/L in our sample, (ii) At low 
stellar M/L the overall relation is driven by age. (iii) For old galax- 
ies there is a trend in dark matter fraction within one R c , from zero 
((M/L) ~ 3) to about 30% dark matter contribution for galaxies 
with the highest M/L in our sample ((M/L) ~ 6). 



5 DISCUSSION 

5.1 A second parameter in the M/L variations? 

In this Section we investigate whether the variations in the M/L 
of the galaxies, apart from a main dependence on the galaxy cr 
(and closely related quantities like mass and luminosity), are also 
dependent on other global observables, in particular on the galaxy 
kinematics, morphology and environment. 

As mentioned in Section Q a number of studies over the 
past twenty years have suggested that early-type galaxies display 
a dichotomy between generally massive, slowly-rotating, metal- 
rich galaxies and rotationally-supported, less massive, metal-poor 
galaxies, which suggests a difference in thei r formation scenar- 
ios (e.g. D avies et alJ ll983|: |Faber et al]|l99 7l). The se differences 
in galaxy structure led Kormendv & Bender] il996f> to propose a 
physically-motivated revision to the classic and stil l broadly used 
(e.g. RC3) visual classification scheme by Hubble ( 1936), based 
on photometry of E and SO galaxies, to take the findings from the 
kinematics into account. The fainter and faster rotating galaxies 
also tend to show steeper luminosity profiles than the brighter ones 
jFaber et alJll997l) and there is indication that they may all con- 
tain disks" TLaueret alj200^) . The SAURON kinematical maps pre- 
sented in Paper III also show two general velocity field morpholo- 
gies: galaxies with a clear sense of rotation, with only small kine- 
matic twists; and galaxies with little or no rotation, or with strong 
kinematic misalignments. The relation between all these different 
galaxy properties is still not fully understood, and it will be ad- 
dressed in a future paper, for the full SAURON representative sam- 
ple. 

Here we limit ourselves to studying the effect on the M / L of a 
single, hopefully representative, parameter: the galaxy kinematics. 
A classic and simple way to quantify the diff erences in gala xy kine- 
matics is by using the anisotropy diagram iBinnev 1978), which 
relates the observed flattening e of a galaxy to the ratio between 
the ordered and random rotation of the s tars (V/o). The formalism 
was recently updated by Binnev 1 2005b) for use with integral-field 
kinematics. The (V/a) — e diagram has often been used to pro- 
vide a separation between fast rotating and slow rotating early-type 
galaxies (e.g. lGeha. Guhathakurta. & van der Marei2003l) . A limi- 
tation of the diagram to characterize the kinematics of a galaxy is 
that it contains no spatial information, so that it cannot distinguish 
between a fast rotating nuclear component and a global galaxy rota- 
tion. To overcome this limitation we introduced a new quantitative 
classification scheme for early-type galaxies which is based on an 
approximate measure of the specific angular momentum of galax- 
ies from integral-field stellar kinematics. In this work we adopt that 
classification to define the galaxies in our sample with significant 
angular momentum per unit mass, which we call the "fast-rotators" 
and the ones with negligible amount of specific angular momentum, 
which we define "slow-rotators" (Table0. The distinction between 
the two classes of galaxies can be easily seen on the kinematic maps 
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of Paper III. More details on the classification will be given in Em- 
sellem et al. (in preparation). 

In the first (top) panel of Fie. ll8l we show the distribution of 
the residuals in the (M/L)-a correlation for the fast-rotators and 
slow-rotators, as defined above. There appears to be a tendency for 
the slow-rotators to have a higher M/L than the fast-rotators, at 
given a c . This difference in the distribution of the residuals between 
the two classes of object s can be q uantified using the Kolmogorov- 
Smirnov (K-S) test (e.g. lPress et"ailll99l Section 14.3). We find 
that the two distributions are significantly different, the probability 
that they are drawn from the same distribution being 2%. The same 
probability is f ound using a com pletely different estimator, the me- 
dian statistics iGott et alj|200ll) . counting the number of galaxies 
above and below the best-fitting relation. 

A similar difference between the M/L of the slow-rotators 
and fast-rotators is seen with respect to the [M/L) pop in Fig. ll7l (or 
Fig. l 161 . in the sense that the slow rotating galaxies tend to have the 
largest M/L at any given (M/ L) pop (or H/3). If one excludes from 
the comparison the five youngest galaxies, and only considers the 
older galaxies for which the (M/L) pop is likely to be reliable, one 
finds that 7/8 of the slow rotating galaxies lie on top or above the 
dashed line (which roughly corresponds to a ~ 30% contribution of 
dark matter within R c ), while only 3/12 of the fast rotating galaxies 
lie above the same line. In this case the probability that this differ- 
ence in M/L is due to a random fluctuation is less than 1%. The 
only slow rotator without evidence for dark matter is NGC 4458, 
which appears peculiar also for being the only low-luminosity and 
low-cr c slow rotator. 

In the second panel of Fig. ll8l we show the dependence of the 
residuals of the (M / L)-a relation on the galaxy luminosity. We 
separated the sample into the galaxies brighter and fainter than the 
K-band luminosity Mk = —24. This ad hoc luminosity value was 
chosen to have the same number of objects in the two classes as 
in the comparison of the top panel between fast and slow rotating 
galaxies. We do not detect any significant difference in the distri- 
bution of the residuals of the bright and faint galaxies. The K-S 
probability that the residuals come from the same distribution is in 
this case 32%. 

This finding that the slow rotating galaxies, which are more 
common among the most massive ones (Fig. I18> . have a higher 
M/L than the faster rotating and generally less massive ones, is 
reminiscent of recent findings, derived using completely differ- 
ent techniques, either from the dynam ics of the planetary nebulae 
iRomanowskv et al l2003t|Napolitan o et alj 2005h or from the anal- 
ysis o f gravitationallv-lensed quasars iFerreras. Saha. & W illiams 
2005), indicating smaller fractions of dark matter among faint el- 
lipticals, compared to more massive ones. However, if our result 
is correct, it would indicate that the dark matter fraction is more 
closely linked to the galaxy dynamics than to the mass alone. 

In the third panel of Fig. l 18l we compare the distribution of the 
residuals in the (M/L)-a relation for the morphologically classi- 
fied E and SO (from RC3). Consistent with previous studies (e.g. 
J0rge nsen et alj[l993) we find no difference in the M/L distribu- 
tion for the two samples. The K-S probability that the two M/L 
distributions are the same is in fact 74%. The lack of any depen- 
dence on galaxy morphology may just reflect the limitations of the 
classic E/SO visual classification, with respect to more physically- 
motivated ones, which als o consider the galaxy kinematics (e.g. 
iKormendv & BendeJl996l) . 

Finally in the fourth (bottom) panel of Fig.^|we test whether 
the M/L shows a different distribution for the subsample of galax- 
ies in clusters and in the field (defined as in Paper II). Hierarchical 
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Figure 18. Residuals in the (M / L)-cr correlation of Fig. [3] The values 
shown are the same in the four panels, but the meaning of the red and blue 
symbols is different in each panel. First Panel: the different galaxies are 
coloured according to their kinematics, using the classification in Table fTI 
Namely the blue diamonds are fast-rotator galaxies, while the red squares 
are the slow-rotator galaxies. Second Panel: the blue diamonds and red 
squares are the galaxies respectively fainter and brighter than Mx = —24. 
Third Panel: here blue diamonds are SO galaxies and red squares are E 
galaxies, as defined in RC3. Fourth Panel: here blue diamonds are field 
galaxies, while red squares are cluster galaxies (as defined in Paper II). The 
K-S probability that the observed distribution of the residuals, for each pair 
of classes of objects, is drawn from the same distribution is 2%, 32%, 74% 
and 10% from the top to the bottom panel respectively. The only signifi- 
cant difference in M/L is thus observed between the fast-rotators and the 
slow-rotators, while we do not detect any significant difference as a function 
of luminosity or morphological classification. There is some very marginal 
evidence for a difference in M/L as a function of environment. 

galaxy formation scenarios jDiaferio et ai]|200ll) predict that the 
galaxies in the field should have a significantly lower M/L than 
the galaxies in clusters, nearly independently of redshift. From our 
sample we find some very marginal difference between the two dis- 
tributions of residuals, the K-S probability that the two distribu- 
tions are the same being 10%. This small difference goes in the 
opposite direction than the model predictions, in the sense that 
the field galaxies have a marginally higher M/L than the cluster 
galaxies. This is consistent with recent findings from high-redshift 
studies, which show a similar evolution for the M/L of cluster 
and field galaxies (e.g. van de Ven. van Dokkum, & Franxl 120031: 
Ivan der Wei et alj|2005t) . Together with the existence of the corre- 
lations we discuss in this paper, this result shows that the physical 
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properties of the individual galaxies play a more important role in 
their formation and evolution than the environment. 

5.2 Origin of the (M/L)-a correlation 

The [M/ L)-a c correlation, which we studied in detail in this paper, 
is part of the remarkable series of tight relations of the a c in early- 
type galaxies and spiral bulges with (i) luminosity, (ii) colours, (iii) 
metallicity and more recently a lso (iv) BH mass and (v) mass of 
the dark halo iFerraresel 120021) . Together with these regularities, 
early-type galaxies are characterized by dichotomies, between the 
most massive giant ellipticals on one side, which are red, metal- 
rich, slowly-rotating, and have shallow photometric profiles, and the 
fainter objects on the other side, which tend to be bluer, metal-poor, 
are dominated by rotation, have cuspy profiles and may c ontain 
disks iDaviesetal. 1983; Kormendv & Bender 1996; Faberetal. 
ll997tlLauer et all2005l) . One way to explain the observed trends is 
with the standard hierarchical scenario, where galaxies are assem- 
bled from smaller building blocks. In this picture the importance 
of gas during the merger process declines with increasing galaxy 
mass, due to the transformation of gas into stars, while galaxies are 
assembled. The role of the gas dissipation in the last merger could 
explain the observed differe nces in the ch aracteristics of early-type 
galaxies I Bender, Burstein, & Faber 1992). However it was realized 
some time ago that an apparently more successful way of explain- 
ing both the regular relations and the dichotomies is to assume some 
form of feedback process lDekel&SilklT 986). which must be able 
to stop the star formation and the chemical evolution by heating and 
expelling the gas from the galaxy potential well (which is related to 
the observed a). 

Until recently the most natural candidate feedback generator 
was the explosion of supemovae. But the last few years of BH 
studies have started to revolutionize this picture. Just five years 
ago the Mbh-? relation first suggested a close connection be- 
tween the formation of the BH and of the galaxy stellar compo- 
nent. Now a number of lines of evidence seem to converge to- 
wards a scenario in which the BH is the main actor in the feed- 
back process which shapes galaxy evolution (e.g. Silk & Rees 1998; 
[ Granato et alj|2004 lBinne^l2005al; lOstriker. & Ciottill2005t Loebl 
l2005t ISprinsel. Pi Matteo. & Hernquistl fcOOSt) . In this picture the 
gas accumulates in the dark matter potential well, feeding the BH, 
forming stars and enriching the medium with metals, until the BH 
grows to the point where its powerful jets heat and expel the gas 
from the galaxy. This scenario naturally explains a close connec tion 
of a with the BH mass I Di Matteo. Springel " & Hernauisl2005h . as 
well as with the asymptotic circular velocity (related to the dark 
halo mass) and the M/L (related to both the halo mass and the 
metal enrichment). 

The most obscure part of this picture is the link between the 
dark and the visible matter, both in the centres and the outskirts of 
early-type galaxies. It is remarkable that, thirty years after the first 
evidence of dark matter in spiral galaxies, a clear evidence of dark 
matter in early-type galaxies is still missing. Although a number of 
isolated studies from stellar dynamical models, X-ray and HI gas 
kinematics, or gravitational lensing indicate, not surprisingly, that 
early-type galaxies are also embedded in massive dark halos, other 
results have challenged this picture, suggesting that some galaxies 
may not contain dark matter after all. In this context our work helps 
to constrain the elusive dark matter contribution in the central re- 
gions of galaxies. The accurate M/ L determinations in galaxy cen- 
tres, compared with M / L derived by kinematical tracers at large 
radii, are also needed to measure dark matter profiles. 



Progress on the investigation of the (M / L)-o relation can 
be made using: (i) larger galaxy samples, to reliably determine 
the intrinsic scatter and the significance of the apparent connec- 
tion between the dynamics and the dark matter content in galaxies; 
(ii) more accurate distance s (e.g. using the SBF method from the 
ACS Virgo Cluster Survey ICote et al.l2 0041. which includes some 
galaxies in our sample); (iii) integral-field observations of samples 
of galaxies in clusters at intermediate distance, where the kinemat- 
ics can still be resolved, but galaxies can be assumed to lie at the 
same distance. 



5.3 Caveats 

The galaxy sample we use in this paper was extracted from the 
SAURON representative sample of E/SO galaxies in Paper II, and 
in addition includes M32 to explore the small mass regime. As de- 
scribed in Section |2~T1 the main selection criteria were the avail- 
ability of HST photometry and SBF distances for the galaxies, but 
we also excluded some obvious bars (a similar selection affects the 
well known Mbh-ct relation). Although we have no indication that 
these selection criteria should bias our results in any specific way, 
they make the sample not necessarily representative any more. The 
fact that we can recover the FP correlations suggests that our small 
sample does represent the E/SO population, but it is clear that in- 
creasing the galaxy modeling sample in a statistically selected way 
would be beneficial for future studies. 

The slow rotating galaxies that we discussed in Section BTI are 
generally thought to be triaxial (e.g. de Zeeuw & Fra nxlll99lh . If 
this is the case our axisymmetric models will only provide a first 
order approximation to the true orbital distribution of the galaxies, 
and the M / L may not be correctly recovered. A significant triaxi- 
ality however may be expected to increase the observed scatter in 
the (A// L)~a relation, which is instead quite small, arguing against 
strong triaxiality of the slow rotating galaxies in our sample. This 
idea is also supported by the fact that these galaxies all appear quite 
round in projection. In any case it is not obvious that axisymmet- 
ric models of a sample of possibly triaxial galaxies, seen at random 
orientations, should produce a positive bias i n the inferred M/L . 
A detailed test using triaxial galaxy models IVerolmeet ai] f2003) 
goes beyond the scope of the present paper. 

In our models we assumed a constant M/L throughout the 
galaxies and in Section l33l we discussed the possible implications 
of this assumption. It is possible for the dark matter profile to vary 
systematically as a function of galaxy mass. As a result our constant 
M/L approximation would not be equally accurate for the galaxies 
in the different mass ranges. This effect could also produce a bias 
on the measured M / L, which is however difficult to quantify, given 
the limited information on the dark matter profiles. 

The result of Section l3~3l shows. however, that it is not easy to 
'hide' dark-matter in the centre of galaxies. For this reason the tight 
( M/ L)-a correlation implies that (i) either dark matter is unimpor- 
tant in galaxy centres or (ii) dark matter is closely related to the lu- 
minous matter. In fact a dark matter halo that dominates the central 
regions more than the one conside red (e.g. the steep er profiles found 
in the numerical simulations of iNavarro. Frenk. & Whitelll99d) . 
would easily erase the correlation derived from the M/L of our 
models, unless the halo properties are closely linked to the stel- 
lar component. A similar conclu sion was reached in a stu dy of the 
tightness of the FP by Borriello, Salu cci. & Danesd J2003). 
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6 CONCLUSIONS 

We measured the M/L in the central regions, for a sample of 25 
early-type galaxies with SAURON integral-field kinematics, using 
both two-integral Jeans models and three-integral Schwarzschild 
models. The galaxy sample was extracted from the SAURON sur- 
vey (Paper II) and spans over a factor of a hundred in mass. To this 
sample we added the galaxy M32 to explore the low-mass range. 
The models are constrained by SAURON kinematics, to about one 
effective (half-light) radius Re, and HST/WFPC2 + MDM ground- 
based photometry in the J-band. We studied the correlations of 
the M/L with global observables of the galaxies and found a se- 
quence of increasing scatter when the M/L is correlated to: (i) the 
luminosity-weighted second moment of the velocity inside the half- 
light radius <j c , (ii) the galaxy mass, (iii) the /-("-band luminosity and 
(iv) the J-band luminosity. These correlations suggest that the M/ L 
depends primarily on galaxy mass, but at a given mass it is a func- 
tion of a, which in turn is related the galaxy compactness. For our 
galaxy sample the (M/L)-a e relation has an observed scatter of 
18% and an inferred intrinsic scatter of ~ 13%. 

We compared the M/L derived from the dynamical models, 
which does not depend on any assumption of spatial or dynami- 
cal homology, with the classic predictions for the (M/L)-pp de- 
rived from the Fundamental Plane, which fully depends on the as- 
sumptions of homology and virial equilibrium. The slope of our 
observed correlations of M/L with luminosity, mass and a, can ac- 
count for ~ 90% of the tilt of the FP. We also compared directly 
the M/L from the models with the virial predictions for our own 
sample, and found a relation (M/L) oc (M/L)^° 8±0 07 , which 
confirms, with smaller uncertainties, the results from the FP com- 
parisons. This also shows that, when the virial parameters R c and 
L are determined in the 'classic' way, using growth curves and the 
<7 C is measured in a large aperture, the virial mass is a reliable 
and unbiased estimator of galaxy mass. This has implications for 
high redshift studies, where the construction of general dynamical 
models is unfeasible. The best-fitting virial relation has the form 
(M/L) v ir = (5.0 ±0.1) x R c ac/(LG), in broad agreement with 
simple theoretical predictions. 

These correlations raise the question whether the observed 
M/L variations are primarily due to differences in the (M/L) pop 
of the stellar population or differences in the dark-matter fraction. 
To test this we compared our dynamical M/L with simple esti- 
mates of the (M /L) pop , derived from the observed line-strengths 
indices of our galaxies. We find that the M/L generally correlates 
with the (M/L) pop , indicating that the variation in the stellar pop- 
ulation, mainly due to a variation in the luminosity-weighted age, is 
an important factor in driving the observed M/ L variation. The rel- 
atively small scatter in the correlation indicates that the IMF of the 
stellar population varies little among different galaxies, consistent 
with results obtained for spiral galaxies by Bell & de J ond 1200 ll) . 
The small scatter also argues against the cus py dark matter dens ity 
profiles predicted by numerical simulations iNavarro et all 1 996 ). 

The accuracy of our M/L determinations allows us to detect 
significant deviations from a one-to-one correlation, which must be 
due to either variations of the IMF among galaxies, or to variations 
in the dark matter fraction within R c . In the latter case the median 
dark matter fraction for our sample is found to be ~ 30% of the total 
mass, broadly consistent wit h earlier findings from dynamics (e.g. 
lGerhardetai1l200ll:l Thomas et al. 2005) and gravitational lensing 
(e.g. lTreu & Koopmansll2004 iRusin & Kochanekll2005l) . We find 
some evidence for the variation in M/L to be related to the dy- 
namics of the galaxies. In particular the slow rotating galaxies in 



our sample, which are more common among the most luminous ob- 
jects, tend to have a higher M/L at given (M/L) pop than the fast 
rotating and generally fainter galaxies. Assuming a constant IMF 
among galaxies, these results would suggest that the slow rotating 
massive galaxies have a higher (~ 30%) dark matter fraction than 
the fast rotating galaxies. We speculate that this difference in M/L 
indicates a connection between the galaxy assembly history and the 
dark halo structure. 
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APPENDIX A: CONSTRAINTS ON INCLINATION FROM 
TWO-INTEGRAL MODELS 

In this Appendix we present the Jeans two-integral models of five 
galaxies in our sample, for which the inclination can be inde- 
pendently determined by non-dynamical arguments. The galaxies 
NGC 524 and NGC 4459 possess very regular dust disks, while 
NGC 4150 has a more disturbed disk, which are all clearly visible 
on th e unsharp-masked HST/WFPC2 images (Fig. 4 of lSarzi et alj 
I2006L Paper V). We fitted an ellipse to the dust morphology, and, 
assuming the disks are intrinsically circular, and in equilibrium in 
the equatorial plane of an oblate spheroid, we inferred an inclina- 
tion of i = 20° ± 5, i = 54° ± 5 and i = 46° ± 2, respectively 
for NGC 524, NGC 4150 and NGC 4459. T he inclination of the gas 
disk in NGC 2974 was carefully studied bv lKrainovic et all 1 2005) 
who derived an inclination of i — 60° ± 3. Finally in the case of 
NGC 4552 there is no evidence of extended patchy dust absorp- 
tion, but this giant elliptical is unlikely to be intrinsically very flat 
(see discussion in Section l3~2l . We expect a broad range of possible 
inclinations, but low values should be excluded. 

In Fig. IA1I we compare the maps of the observed projected 
second moment of the velocity fi' 2 = VV 2 + a 2 , as derived from 
the SAURON stellar kinematics, with the maps of the unique second 
velocity moment associated to the MGE density distribution, com- 
puted by solving the Jeans equations, at every inclination. To help 
the visual comparison between data and models, the observed fi' 2 
was bi-symmetrized by replacing the value in each bin with the av- 
erage of the corresponding values in the four quadrants (computed 
using linear interpolation). Each plot is scaled to the best-fitting 
M/ L using equation J2j. This same kind of fit was performed for 
all the 25 galaxies in the sample, to determine the inclination at 
which the M/L of the Jeans and Schwarzschild models was mea- 
sured. The Jeans models always provide a clear best fitting value for 
the inclination, which can also be easily recognized by eye. Start- 
ing from an edge-on model (i = 90°), the effect of lowering the 
inclination (while keeping the M/L at the best fitting value) can 
be mainly described as a decrease of the second moments along 
the minor axis, accompanied by a corresponding increase along the 
major axis. For the five galaxies considered in this Appendix the 
best-fitting inclination derived from the Jeans models agrees with 
the expected values. 



APPENDIX B: MGE PARAMETERS OF THE SAMPLE 

In this Appendix (Tables B1-B2) we present the (distance inde- 
pendent) numerical values of the analytically-deconvolved MGE 
parametrization of the surface brightness of the galaxies in our 
sample. 7 The constant-PA models were obtained by fitting the 
HST/WFPC2/F814W photometry at small radii and the ground- 
based MDM photometry (in the same band) at large radii. They 
provide an accurate description of the /-band surface brightness of 

7 The values fo r the MGE parameterization of NGC 2974 were given in 
Krainovic et al. 1 2005) 
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Figure Al. Top Panels: Bi-symmetrized and linearly interpolated map of the observed second velocity moment p' 2 = \JV A + a 2 , extracted from the SAURON 
stellar kinematics of the galaxies NGC 524, NGC 2974, NGC 4150, NGC 4459 and NGC 4552. The actual Voronoi bins in which the kinematics was measured 
are shown in Paper III. The black central disk shows the region that was excluded from the following fits. Subsequent Panels: two-integral Jeans modeling of 
the unique projected second moments associated to the MGE mass density distribution, for different inclinations. Each plot is scaled to the best-fitting M/L 
using equation |2). For every galaxy the lowest inclination corresponds to an intrinsic axial ratio q = 0.1 of the flattest Gaussian component in the MGE 
model. Subsequent panels are equally spaced in q. The inclination of each model, and the corresponding best-fitting M/ L is written above each plot. The 
inclination independently inferred from the geometry of the dust or from the gas kinematics is i m 20° , 60° , 54° , 46° , for NGC 524, NGC 2974, NGC 4150 
and NGC 4459 respectively. No strong constraint exists for the inclination of NGC 4552, but this giant elliptical is not expected to be flat (see Section l3*2*l . so 
it should be best fitted for i PS 90°. For all the five galaxies the best fitting inclination derived from the Jeans models agrees with the expected value, as can 
also be seen by eye on these maps. 



the galaxies from R 
in each galaxy (usually corresponding to 5-10 R c ). The nucleus of 
NGC 4150 is heavily obscured by dust. For this galaxy we fitted the 
nuclear MGE model to an HST/NICMOS/F 1 60W (H-band) image, 
which was rescaled to the WFPC2/F814W (J-band) one for cali- 
bration, neglecting possible I — H color gradients. In other galax- 
ies dust and bright foreground stars were simply excluded from the 
MGE fit. The matching of the different images and the calibration 



adopted for the MGE models is described in Section [3~TI The de- 
convolved MGE surface brightness E is defined as follows: 



N 

£ 

2=1 



Ij exp 



2a] 



i2 . y 

X + —p 



(Bl) 



where the model is composed by TV Gaussian components of dis- 
persion <jj, axial ratio q'j and peak intensity Jj. The coordinates 
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Table Bl. MGE parameters for the deconvolved /-band surface brightness. 
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(x' , y') are measured on the sky plane, with the x'-axis correspond- 
ing to the galaxy major axis. The total luminosity of each Gaussian 
component is given by Lj = 2irl'jO?jqj (see lCap pellari 20o3. for 
details). 
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Table B2. MGE parameters for the deconvolved /-band surface brightness. 



3 


lOg 1 

° J 


logo-j 


''j 


1 „„. 77 

log/j 


log <Tj 


q j 


log 1 


log (Tj 


q j 


log 1 

6 J 


log (Tj 


„/ 
q j 




( L 0,-T P C ~ 


z ) (arcsec) 




(L Q j pc 


*) (arcsec) 




( L 0,/ P c " 


z ) (arcsec) 




( L 0,/ P C " 


z ) (arcsec) 








NGC 4459 






NGC 4473 






NGC 4486 






NGC 4526 




1 


5.592 


-1.737 


0.820 


4.346 


-1.762 


0.550 


6.541 


-1.762 


0.950 


5.463 


-1.230 


0.900 


2 


5.463 


-1.507 


0.820 


3.900 


-1.137 


0.545 


4.738 


-0.971 


0.955 


4.610 


-0.441 


0.900 


3 


4.947 


-0.928 


0.817 


4.277 


-0.623 


0.555 


3.513 


-0.415 


0.944 


4.223 


-0.067 


0.900 


4 


4.560 


-0.513 


0.929 


4.207 


-0.315 


0.528 


3.288 


0.007 


0.990 


4.049 


0.470 


0.474 


5 


4.504 


-0.197 


0.949 


4.134 


-0.021 


0.491 


3.375 


0.520 


0.990 


3.748 


0.629 


0.651 


6 


3.731 


0.115 


0.933 


3.755 


0.212 


0.750 


3.468 


0.834 


0.990 


3.379 


0.969 


0.678 


7 


3.706 


0.370 


0.822 


3.883 


0.306 


0.400 


3.209 


1.114 


0.967 


3.003 


1.329 


0.635 


y 

o 


3.ZZ3 


U. JyO 


n Qon 

VJ.yyVJ 


3.UJ4 


Lf. J 34 


U.OoO 


Z.oOZ 


L.Joy 


f) O/IQ 
xj.y^+y 


1 1 OS 


1 .Do3 


u.vuu 


9 


3.362 


0.765 


0.879 


3.529 


0.583 


0.400 


2.460 


1.731 


0.886 


2.253 


1.718 


0.273 


10 


2.968 


1.021 


0.830 


3.372 


0.833 


0.537 


1.671 


1.906 


0.990 


2.169 


1.938 


0.262 


1 1 


2.631 


1.303 


0.795 


3.078 


0.849 


0.750 


1 .67 1 


2.166 


0.700 


1.594 


2.048 


0.438 


12 


2.109 


1.638 


0.750 


2.922 


1.153 


0.541 


1.153 


2.515 


0.700 


0.558 


2.264 


0.900 


13 


1.632 


1.891 


0.990 


2.483 


1.421 


0.532 














14 








2.149 


1.662 


0.572 














15 








1.639 


1.989 


0.654 


















NGC 4450 






NGC 4552 






NGC 4621 






NGC 4660 




1 


5.384 


-1.596 


0.600 


5.064 


-1.762 


0.990 


5.797 


-1.397 


0.800 


5.984 


-1.620 


0.900 


2 


4.724 


-1.038 


0.600 


4.282 


-0.396 


0.800 


5.323 


-1.030 


0.798 


5.138 


-1.071 


0.900 


3 


4.181 


-0.643 


0.600 


4.325 


-0.339 


0.990 


4.956 


-0.701 


0.773 


4.853 


-0.763 


0.533 


4 


4.067 


-0.501 


0.372 


4.209 


-0.057 


0.957 


4.557 


-0.415 


0.950 


4.709 


-0.541 


0.671 


5 


3.891 


-0.234 


0.600 


4.153 


0.223 


0.933 


4.200 


-0.159 


0.321 


4.336 


-0.209 


0.546 


6 


3.416 


0.208 


0.554 


3.784 


0.594 


0.948 


4.138 


-0.062 


0.893 


3.997 


-0.182 


0.900 


7 


3.046 


0.647 


0.600 


3.416 


0.860 


0.975 


3.838 


0.241 


0.819 


4.187 


0.154 


0.450 


8 


3.218 


0.815 


0.289 


2.711 


1.187 


0.887 


3.900 


0.295 


0.295 


3.946 


0.346 


0.590 


9 


2.847 


1.053 


0.385 


2.608 


1.389 


0.928 


3.607 


0.589 


0.743 


3.319 


0.590 


0.730 


i n 
1 U 


Z. 1 /U 


1 1 1 "7 
1 . J 1 / 


U.ZZo 




1. IZj 


n oo/i 
U.8o4 


1 1 1 "7 

3.31 / 


U.DJO 


U.3UZ 




u.oy / 


U.VUU 


ii 


2.487 


1.381 


0.320 


1.540 


2.136 


0.990 


3.319 


0.886 


0.755 


3.163 


0.981 


0.450 


12 


1.866 


1.480 


0.394 


— 


— 





2.538 


1.094 


0.333 


2.152 


1.148 


0.900 


13 


1.315 


1.686 


0.183 








2.885 


1.277 


0.635 


2.417 


1.265 


0.450 


14 


1.029 


1.867 


0.600 








2.461 


1.544 


0.628 


1.930 


1.451 


0.779 


15 


1.007 


1.867 


0.153 








1.639 


1.812 


0.926 


1.243 


1.806 


0.878 


16 














1.706 


1.849 


0.472 








17 














1.379 


2.190 


0.950 












NGC 5813 






NGC 5845 






NGC 5846 






NGC 7457 




1 


3.645 


-0.457 


0.980 


5.559 


-1.762 


0.800 


3.647 


-1.762 


0.990 


6.308 


-1.762 


0.900 


2 


3.913 


-0.262 


0.980 


5.229 


-1.271 


0.800 


3.794 


0.167 


0.990 


4.778 


-1.059 


0.900 


3 


3.848 


-0.026 


0.900 


4.688 


-1.049 


0.462 


3.125 


0.449 


0.990 


4.056 


-0.680 


0.900 


4 


3.755 


0.254 


0.885 


4.400 


-0.567 


0.800 


3.334 


0.725 


0.990 


3.702 


-0.424 


0.900 


j 


3 401 


522 


924 


4 197 


-0 184 


768 


2 772 


1 063 


v.yyjy 


J.OJJ 


-0 132 


696 


6 


2.962 


0.810 


0.909 


4.232 


0.015 


0.263 


2.625 


1.289 


0.974 


3.451 


0.172 


0.703 


7 


2.516 


1.024 


0.882 


3.862 


0.063 


0.708 


1.946 


1.526 


0.850 


3.174 


0.495 


0.613 


8 


2.264 


1.336 


0.784 


3.710 


0.237 


0.668 


1.984 


1.777 


0.924 


2.541 


0.812 


0.528 


9 


2.180 


1.623 


0.721 


3.629 


0.459 


0.800 


1.431 


2.197 


0.850 


2.535 


0.939 


0.655 


10 


0.706 


2.063 


0.980 


2.632 


0.664 


0.538 








2.520 


1.256 


0.555 


11 


1.467 


2.063 


0.720 


2.909 


0.820 


0.564 








2.236 


1.628 


0.517 


12 








2.127 


1.206 


0.708 








1.369 


1.958 


0.774 



© 2006 RAS, MNRAS OQO.HTIBI 



